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Deep Survey (SXDS) Field. We obtain a photometric sample of 858 LAE candi- 
dates based on deep Subaru/Suprime-Cam imaging data, and a spectroscopic 
sample of 84 confirmed LAEs from Subaru/FOCAS and VLT/VIMOS spec- 
troscopy in a survey volume of ~ 10^ Mpc^ with a limiting Lya luminosity 
of ~ 3 X 10^^ erg s~^. We derive the LFs of Lya and UV-continuum (~ 1500 A) 
for each redshift, taking into account the statistical error and the field-to- field 
variation. We find that the apparent Lja LP shows no significant evolution be- 
tween z = 3.1 and 5.7 within factors of 1.8 and 2.7 in L* and 0*, respectively. 
On the other hand, the UV LP of LAEs increases from 2; = 3.1 to 5.7, indicat- 
ing that galaxies with Lya emission are more common at earlier epochs. We 
identify six LAEs with AGN activities from our spectra combined with VLA, 
Spitzer, and XMM-Newton data. Among the photometrically selected LAEs at 
z = 3.1 and 3.7, only ~ 1 % show AGN activities, while the brightest LAEs with 
logL(Lya) > 43.4 — 43.6 erg s"^ appear to always host AGNs. Our LAEs are 
bluer in UV-continuum color than dropout galaxies, suggesting lower extinction 
and/or younger stellar populations. Our stacking analyses provide upper limits 
to the radio luminosity and the fueii/ fhya line fraction, and constrain the hidden 
star formation (+ low-luminosity AGN) and the primordial population in LAEs. 

Subject headings: galaxies: formation — galaxies: high-redshift — cosmology: 
observations 



1. Introduction 



Lya emitters (LAEs) provide clues to evolution and formation of galaxies. Most of LAEs 
are known as star-forming galaxies at 2 ~ 2 — 7 with a faint-UV continuum, but a prominent 
Lya emission line, which are pro duced by star-forming activities with a t ypical star-formation 



rate of ~ 1 - lOM ^^yr"^ (e.g., ICowie fc Hul Il998l : lOuchi et all l2003l : lOawiser et all I2OO6 



Pirzkal et al.ll2006l ). LAEs are believed to have a negligible fraction of AGN activities, s ince 



an X-ray detection ha s been reported for only one LAE at 2; ~ 3 (IGawiser et al.ll2006l . see 



also IWang et al.l 120041 ). The strong Lyg emission and a blue UV continuum color imply a 



young metal-poor star-forming galaxies (IMalhotra &: Rhoadd 120021 : iPinkelstein et al.l 120071 ). 



Recent infrared observations and population synthesis models indicate that typical stellar 
mass of LAEs is as small as ~ 10^ — lO^Mf;^ and that the stellar age is a s young as < 10 



Myr ( 


Pirzka 


et al. 


Nilsson et al. 


2007) 



see also I Chary et al.l l2005l : iGawiser et al.l l2006l : iLai et al.l 12006 



20071). LAEs show a compact UV morphology with a size smaller than 1 kpc 



( jPascarelle. Windhorst. Keel, fc Odewahnlll996l : iPirzkal et al.ll2006l ). The fraction of LAEs 
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(or galaxies with a Lyg emission line) in creases at faint UV magnitudes (jOuchi et al.l 12003 
Ando et al.ll2006l : IShimasaku et al.ll2006l ) . These pieces of evidence suggest that L AE samples 
include a young low-mass population at high-z. 

More interestingly, some of LAEs have a Lya emission line with a large rest-frame 
equivalent width (EW) of EWo > 240A that cannot be explain e d by a star-formation with 
a Salpeter IMF jMalhotra fc Rhoadsl [2OO2I : IPawson et aPbooi IShimasaku et alibooeh . It 
is reported t hat about 10-40% of spectroscopicall y-identified LAEs have an EWq of > 240 A 
at z = 4.5 (jPawson et al.l 120041 ) and z = 5.7 (IShimasaku et al.l l2006l ). The population 
synthesis models indicate that such large EW obj ects would have a top heavy IMF, a very 
poor metallicity and /or a ve ry young age< 10^ yr (jCharlot fc Falllll993l : iMalhotra fc Rhoads 
2002; Schaerer 2003, see also Kudritzki et al.|[2000 ). Large- EW objects often sho w a spatially- 



exten ded Lya envelope whose projected size can be as large as > 100 kpc (ISteidel et al. 
2OOOI ). These objects are referred to as Lya blobs. Lya blobs are o riginally discovered 
in an overdensity re gion of Lyman break ga laxies (LBGs) at 2; ~ 3 (ISteidel et al.l 12000 



Matsuda et al.ll2004j . see alsolDev et al. |2005). However, r ecent deep surveys also find Lya 
blobs in blank fields flNilsson et al.l 1200 6: Sait o et al.ll2006l ). Some of these Lyg blobs show 



evidence of AGN activities (jPey et al.i i2005 ). while others are not (IMatsuda et al.l 12004 



Saito et al.ll2006l ). Theoretical studies predict that th ese large-EW obje cts with an extended 



Lya profile a re candidates o f cooling clouds (e.g. lYang et al.l l2006l ) and population III 
galaxies (e.g. ISchaered l2003l ). which are at the very beginning stage of galaxy formation. 
Thus, LAEs are keys to understanding the early stage of galaxy formation with less-massive 
(or dwarf) population some of which include young galaxies even at the very beginning 
stage of radiating clouds and population III. T he studies of LAEs complement the recent 



work of high-z ma ssive and/or old galaxies (e.g. Ivan Dokkum et al.ll2006l : iKriek et al.ll2006 
Daddi et al.lboOTh . 



LAEs are usually ide ntified by a red s hifted Lya ernissiori line that falls in a passband 



2004; van Breukelen et al. 2005; Martin et al 



of nar row-band filter ( e.g. IHu et al.ll2002l : iKodaira et al.ll2003l ). tuned at z ~ 3 up to 2; ~ 9 



Ive et all l20od : Iwillis et al.l I2OO6I). Recent blank -field shtless, multi-si l t, and IFU 



(e.g. 

spectroscopy also searched for LAEs (iMartin &: Sawickil l2004l : iTran et al.l l2004l : iKurk et al. 



20061 ). Blind spectroscopy in cluster regions 



yielded candidates of gravitationally lensed LAEs that enable us to study intrinsically faint 
LAEs JeUIs et allboOll : ISantos et al.l 12004 : Istark et allboO?! ). The numbers of p hotomet 
rically and spectroscopi c ally identified LAEs are riow over several hundreds (e.g. IHu et al. 



2002 




2006 





Ouchi et al. 



2003; 



Kashikawa et al. 



Mal 



2006 



lotra fc Rhoads 



2004 



Venemans et al. 



Taniguchi et al.l l2005l : IShimasaku et al. 



20071 ). an d the numbers are i ncreas ing by re 



cent and on- going w i de- field projects such as WFILAS (jWestra et a. 



(iGawiser et al.ll2006l : iGronwall et al.ll2007l ). Hawaii (see lHu fc Cowie 



20051, 12OO6D, MUSYC 



2006f ). and SSA22 sur- 
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veys (Nakamura et al. in preparation). 

In spite of increasing observational data of LAEs, evolution of LAEs is not clearly under- 



at 2; = 3 — 6 within the errors of measurements, (e.g. 


Rhoads & Mai 


lotra 


2001; 


Ouchi et al. 


2003; 


Hu et a 




2004 


: Yamada et al. 


2005; 


van '. 


3reukelen et al. 


2005 


Shimasaku et al. 


2006; 


Tapk 


en et al. 


200d: ] 


Uuravama et al. 


2007 


, c.f. 


Maier et al. 


2003). Although these measure- 



ments qualitatively agree, the constraints on the LP evolution are weak, due to the small 
statistics and f ield variations. Indeed, the presence of l arge-scale structures of LAEs has 
been reported (IShimasaku et al.ll2003l ; lOuchi et al.ll2005a) at these redshift s , which produc e 
a significant source gradient even within a 0.2 deg^ field ( lOuchi et al.ll2003l ; IHu et al.l 120041 ). 
Note that the narrow-band LAE studies observe a thin slice [Az = 0.1) of the Universe 
whose survey volume is about 1/10 of th at of dropout or LBG surveys {Az = 1) at a given 
search area (see e.g. lOuchi et al.ll2004bl ). We need to obtain statistical properties of LAEs 
including clustering with a good accuracy. Moreover, we should investigate rare populations 
of LAEs with a large EW or AGN, which are essential for understanding galaxy formation 
and its relation to the formation of massive black holes. 

Here we carry out a systematic wide-field narrow-band survey for LAEs at z = 3.1 — 6.6 
in the wide-field (1 deg ^) Subaru/XMM-Newton Deep Field (SXDS) that is covered by 
deep Subaru optical images (^'[5cr] = 27.0; Furusawa in p reparation) and the n ear-infrared 
images of UKIDSS /Ultra D eep Survey (UKIDSS/UDS; [Lawrence et al.l l2006l ) as well as 
VLA (Simpson et al.l 2006a), XMM-Newton (Ueda et al. in preparation), and Spitzer 
( iLonsdale et al.ll2003l ). Except for the data of z = 6.6 LAEs, we have completed narrow- 
band imaging for z = 3.1, 3.7 and 5.7 LAEs. In this paper, we show properties of LAEs 
at z = 3.1 — 5.7. Note that the first results from the z = 5.7 data have been presented in 



Ouchi et al.l (j2005al ). Our systematic survey will not only put strong constraints on prop- 
erties of LAEs at 2 = 3.1 — 5.7, but also give the useful basel ines to be compared with 



higher redshift LAEs for studying reionization of the Universe ( iMalhotra &: Rhoadd 12004 
Kashikawa et al.ll2006l ). 



We show photometric and spectroscopic data, and make the spectroscopic samples in 
Section [2J We define the photometric samples of LAEs in Section [31 The AGNs in our 
LAE samples are presented in Section HI We derive LFs of LAEs and investigate evolution 
of the LF in Section [51 In Section [6l we address the properties of LAEs, such as EW, 
UV-continuum slope, star-formation rate, and AGN activity. We discuss our results and 
summarize them in Sections [7] and [HI respectively. We will present clustering and stellar 
population of these LAEs in two companion papers of Ouchi et al. (in preparation) and Ono 
et al. (in preparation). 
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^Throughout this paper, magnitudes are in the AB system (jOkdll974j : iFukugita. Shimasaku. fc Ichikawa 

19951 ). EW is presented in rest frame (i.e. EWq). The values for the cosmological parameters 



adopted in this paper are: {h, Q^, ^a, n) = (0.7, 0.3, 0.7, 1.0). 



Observations and Data Reduction 



2.1. Imaging Data 



W e carried out na rrow-band imaging with Subaru/Suprime-Cam (IMiyazaki et al.l 12002 



see also llye et aLll2004l ) in 7 nights during 2003 September 28-30, October 22-24, and October 
26. We summarize details of these observations as well as image qualities in Table [H We 
obtained narrow-band images with three bands, NB503 (A^ = 5029A, AA = 74A), NB570 
(Ac = 5703A, AA = 69A), and NB816 (A^ = 8150A, AA = 120A). We show the response 
curves of our narrow-band and broad-band filters in Figure [H These response curves include 
atmospheric absorption, quantum efficiency, and transmittance of optical elements of the 
telescope and instrument. We have chosen these three narrow bands to identify LAEs at 
z ~ 3.1, 3.7, and 5.7. Note that there are no strong OH sky lines within the passbands of these 
narrow-band filters (Figured]). Since Suprime-Cam has a field-of-view (FoV) of 0.255 deg^, 
the 1 deg^ area of SXDS is covered with 5 pointings of Suprime-Cam. These 5 pointings 
are referred to as SXDS-C (Center), SXDS-N (North), SXDS-S (South), SXDS-E (East), 
and SXDS-W (West) whose central coordinates are the same as the archival broad-band 
images of the SXDS project (see Table 1 of Furusawa et al. in preparation). The on-source 
exposure times of A^i?503, NB570, and NB816 filters are typically 1.2 — 1.5, 1.2 — 1.5, and 
4.0 — 5.7 hours per pointing, respectively. The observational condition was photometric and 
clear throughout the first 6 nights, and especially good in three nights for NB816 data in 
2003 September 28-30. The seeing size was typically 0".4 - 0".6 and 0".5 - 1".0 for 2003 
September 28-30, and October 22-24 and 26, respectively. In addition to these narrow-band 
data, we use archival data of very deep broad-band {B, V, R, i' and z') images of the SXDS 
project. The narro^ y-band data are reduced with the S uprime-Cam Deep field REDuction 
package (SDFRED: lYagi et al.ll2002l : lOuchi et al.ll2004al ). With the standard parameter sets 
of SDFRED, we perform bias subtraction, fiat-fielding, distortion+atmospheric-dispersion 
corrections, sky subtraction, image alignments, and stacking. Before stacking, we mask out 
bad data areas such as dead pixels and satellite trails. Cosmic rays are removed in the process 
of stacking with the rejected-mean algorithm. The 5a sky noise of the reduced narrow-band 
images are iVS503_~ 25.3, NB570 ~ 24.9, and NB816 ~ 26.0 in a 2".0-diameter circular 
aperture (see 

reduction are ~ 0".8 in FWHM. 



Ouchi et al.l l2005al for the A^S816 data). Typical PSF sizes of images after 
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We do not use, either contaminated areas with halos of bright stars and CCD blooming 
or low signal-to-noise (S/N) regions located around the edge of the FoV, which are caused by 
dithering. After we reject these bad areas, the effective total areas are 0.983, 0.965, and 1.034 
deg^ for NB503, NB570, and NB816 images, respectively. These effective areas correspond 



to the survey volumes of 7.0 x 10^ Mpc^ 



3.1), 6.1 X 10^ Mpc=^ 



3.7), and 9.2 x 10^ 



Mpc^ 



5.7), if we assume a simple top-hat selection function of LAEs whose redshift 



distribution is defined by the FWHM of narrow-band filters. 

During the observations, we took images of spectrophotometric standard stars of GD50 
and GD71 in NB503 and NBS16 bands and GD248, G93-48, and GD50 in NB570 band 



( lOkdll990l : iBohlin et al.lll995l ). These standard stars were observed 1 — 4 time(s), when the 
sky was thought to be photometric. We calculate photometric zero-points from photometry 
of standard stars. We check these photometric zero points with the narrow- and broad- 
band images by comparing with colors of stellar objects in our field and 175 Galactic stars 
calculated from spectra given in Gunn fc Strykei (Il983). We find that the colors of stellar 
objects in our data are consistent with those of iGunn fc Strykerl 's stars within 0.01 — 0.05 
magnitude. The photometric-zero points thus obtained are regarded as more accurate than 
0.05 mag in the entire 1 deg^ field. 

Our narrow-band images are registered to match the coordinates of SXDS broad-band 
images based on hundreds of stellar objects commonly detected in both narrow- and broad- 
band images. The astrometry of our objects is the same as those of SXDS version 1.0 catalog 
(Furusawa et al. in preparation). The errors in the relative positions of objects are ~ 0".04 
in r.m.s. The r.m.s accuracy of the absolute positions is estimated in Furusawa et al. (in 
preparation) to be ~ 0".2. After the registration, we homogenize the PSF sizes of broad and 
narrow-band images, referring to these stellar objects. The PSF sizes of narrow-band images 
match to those of broad-band images with an accuracy of 5FWHM ~ 0".01. 



2.2. Photometric Catalogs 



Source detection and photometry are performed using SExtractor (IBertin fc Arnouts 



19961 ). We measure both MAG_AUTO of SExtractor and 2".0-diameter aperture magni- 
tudes. We adopt MAG_AUTO as total magnitudes, while we use a 2".0-diameter aperture 
magnitude to measure colors of objects in order to obtain colors of faint objects with a good 
signal-to-noise ratio. We make A^-B503-, NB570-, and A^i?816-detection catalogs, and limit 
to NB503 < 25.3, NB570 < 24.7, and NB816 < 26.0, respectively, that correspond to 
about 5a detection limits on a 2".0-diameter aperture. Our A^i?503-, NB570-, and A^i?816- 
detection catalogs contain 98,907, 64,362, and 278,458, respectively. We correct the magni- 
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tudes of objects for Galactic extinction of E{B — V) = 0.020 (ISchlegel. Finkbeiner. fc Davis 
1998h . 



2.3. Spectroscopic Data 



We carried out spectroscopic follow-up o 
Object Camera and Spectrograph (F OCAS; 



Dservations for our LAE candidates with Faint 



Multi-Object Spectrograph (VIMOS; iLeFevre et al. 



Kashi 



cawa et al.l |2002| ) of Subaru and Visible 



20031) of VLT. 



The FOCAS observations were carried out in the MOS mode on 2003 December 20, 23, 
and 25, 2004 October 17, November 9, 11, and 2005 November 2-3. The sky was clear during 
these observations, except on 2005 November 3. Since these observations were conducted 
under the SXDS project and the Subaru open use programs of Akiyama et al., Yamada et 
al., and Sekiguchi et al., the slits of our objects shared the eight MOS masks with those of 
the other targets. The on-source exposure time ranges typically from 7200 to 10800 seconds. 
These spectroscopic observations were made with SY47 order-cut filter and 300 lines mm~^ 
grating having the blaze wavelengths of 5500A (300B) with a slit width of 0".8. The spectral 
range and resolution are A = 4900 - 9400A and A/ A A ~ 500 - 1000, respectively. 

The VI MOS spectroscopy was conducted in the programs of Simpson et al. (in prepa- 
ration) and ISaito et al.l (120071 ). The Simpson et al.'s observations were made between UT 
2004 Dec 17 and UT 2006 Jan 2. The medium resolution grating was used with the GG475 
order sorting filter, which provides a spectral resolution A/ A A ~ 580 over the wavelength 
range 4800 A-1 /im. Each mask was observed with 1 x 2700 + 2 x 1350 second exposures, 
often on separate nights. Data reduction is described in detail in Simpson et al. (in prepa- 
ration) but broadly followed the standard pipeline method. The Saito et al.'s observations 
were carried out in 2004 November 6 to 9. The HR-Orange grism and the GG435 order-cut 
filter were used for MOS masks with a slit width of 1".0. The spectral resolution of these 
observations is i? ~ 2160A, which is > 2 — 4 times higher than the other observations of our 
spectroscopy. The effective on-source i ntegration t i me w as 4-7 hours. Data reduction and 
details of observations are described in lSaito et al.l (120071 ). 



2.4. Spectroscopic Catalogs and Samples 

We took spectra for objects with colors similar to LAEs which show a narrow-band excess 
and a continuum break in blue band. We observed 128 and 29 objects at z = 3.1/3.7 and 5.7, 
respectively. Since we share the MOS masks with targets from the other projects, we had 
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tight constraints on our target selection for spectroscopy. Thus, we include unreliable LAE 
candidates with colors similar to foreground objects and very faint narrow-band magnitudes, 
some of which potentially have a bright Lya emission line redshifted to the edge of passband 
of our narrow band. As a result, we have identified line emitters from 60% of our targets. 
Although the success rate of identification is not high by the constraints of our spectroscopy, 
this target selection allows us to obtain spectra for objects in wide color and magnitude 
ranges. 

In addition to these spectra for our LAEs, we use the SXDS version 1.0 spectroscopic 
data taken with Subaru/FOCAS and AA0/2dF. The combination of these and our data 
provide 3,233 spectroscopic redshifts in our field (Akiyama et al. in preparation). We refer 
this catalog to estimate the contamination rates of our LAE sample in Section 13.21 Figure [2] 
presents the redshift distribution of our sources from the combined catalogs. 

We have investigated our spectra with an emission line at the passbands of our narrow 
bands. We carefully check the possibilities of low- 2; [Oil], [Olll], and Ha emitters. For 
z = 3.1 LAE candidates with an emission at ~ 5000 A, we look at the spectrum at the 
wavelengths of an Ha line of a 2; = 0.004 [Olll] emitter and Ha/[Olll] lines of a 2; = 0.3 
[On] emitter. Similarly, the possibilities of a z = 0.1 [Olll] and a z = 0.5 [Oil] emitters are 
examined for z = 3.7 LAE candidates. For z = 5.7 LAE candidates, our spectra allow us to 
identify a. z = 0.2 Ha emitter and a z = 0.6 [Olll] emitter with corresponding [Oll]/[Olll] 
and [On] lines, respectively. However, no nebular emission lines enter in the wavelengths of 
our optical spectra for a z = 1.2 [Oil] emitter. Because our spectral resolution is not high 
enough to clearly identify an [Oil] doublet, we examine the possibility of an [Oil] emitter 
by the detection of blue continuum in S or bands whose 2a upper limits reach very deep 
magnitudes oi B = 28.7 and V = 27.7. In this way, we discriminate low-z emitters from 
high-2; LAEs. We also visually inspect images and spectra, and remove spurious objects. 
By these analyses, we have identified 41, 26, and 17 LAEs at z = 3.1, 3.7, 5.7. We refer 
to these 84 LAEs as the spectroscopic samples of LAEs. In Figures [3ll5] we show spectra of 
z = 3.1 — 5.7 LAEs, together with snapshots of broad- and narrow-band images. Table [2] 
summarizes the properties of LAEs with a spectrum. Figure [6] plots the redshift distribution 
of these spectroscopically-identified LAEs at z = 3.1, 3.7, and 5.7. 
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3. Photometric Samples of Lya Emitters at z = 3.1 — 5.7 

3.1. Definitions oi z = 3.1, 3.7, and 5.7 Lya Emitters 

We plot color- magnitude diagrams in Figure [7] for objects in our photometric and spec- 
troscopic catalogs. Each panel of Figure [7] shows narrow-band excess color and narrow-band 
magnitude for NB503, NB570, and NB816. Figure [8] presents two-color diagrams based on 
the NB503-, NB570-, and iVi?816-detection catalogs, together with colors of model galaxies 
and Galactic stars. Colors of model objects indicate that LAEs can be isolated from low-z 
galaxies and Galactic stars by their narrow-band excess of Lya emission and red continuum 
colors. 

We compare colors of galaxies in our LAE samples with those of the 3,233 spectroscopically- 
identified objects which include LAEs and foreground/background interlopers. As expected, 
spectroscopically-identified LAEs are located in the upper-right part of the two-color dia- 
grams. 

Based on these color diagrams, we define the selection criteria of three LAE samples: 

V - NB503 > 1.2 & [{V <V2^ &^ B -V > 0.5) or (V" > & S - > 0.5)] for z = 3.1 LAEs, (1) 

V - NB570 > 1.3 & [(y <V2a ^ B-V> 0.7) or (V >V2a ^ B - ¥2^ > 0.7)] for 3.7 LAEs, (2) 
i' - NBSm >1.2 k B > B2a kV > ¥2^ & [{R <R2a k R-i' > 1.0) or (i? > i?2^)] for z^ 5.7 LAEs, (3) 

where i?2CT, V2a, and are 2a limiting magnitudes of B, V, and R images, i.e. B = 28.7, 
V = 27.7, and R = 28.0, respectively. With these criteria, we include the objects with 
no detectable continuum in V (for z = 3.1 and 3.7) and i' (for z = 5.7). Thus, these 
photometric samples are complete to the narrow-band magnitude limits, which include all the 
Lya emitting objects at these redshifts either with or without a detectable-UV continuum. 
The narrow-band excess colors, V-NB503 > 1.2, V-NB570 > 1.3, and i'-NB81Q > 1.2, 
isolate an emission line object at each redshift with a fiat continuum {f,y =const) whose rest- 
frame EW is greater than ~ 45A. Note that the actual EW limits for our LAE samples are 
different from ~ 45A, since LAEs do not have a flat continuum, but a continuum break, i.e. 
Gunn-Peterson trough, between blue and red sides of Lya emission. Thus, the narrow- and 
broad-band colors do not provide a clear limit of EWs. We estimate the limits of rest-frame 
EWs to be approximately ~ 64A, ~ 44A, ~ 27A for our z = 3.1, 3.7, and 5.7 LAE samples 
with models of our Monte-Carlo simulations described in Section 15.11 

We apply these selection criteria to our photometric catalogs. We use A^-B503-, NB570-, 
and iVi?816-detection catalogs and flnd 356, 101, and 401 LAEs at 2; = 3.1, 3.7, and 5.7, 
respectively. In this way, we obtain LAE samples composed of 858 photometrically-identifled 
LAEs in total. We refer to these 858 LAEs as the photometric samples of LAEs, which are 
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summarized in Table O Note that these photometric samples are purely composed of objects 
satisfying the selection criteria of eq. ([1]), ([2]), or ([3]). Our photometric samples do not include 
spectroscopically-identified LAEs escaping from the photometric criteria. 

Figures [QlfTTl show surface densities of LAEs. The red circles are the average surface 
densities. The black points with 5 different symbols indicate the surface densities in 5 sub- 
fields (~ 0.2 deg^) of Suprime-Cam, i.e., SXDS-C, -N, -S, -E, -W. The detection completeness 
is corrected with the simulations described in Section 13. 2[ For these completeness-corrected 
points, the difference of surface densities among the five ~ 0.2 deg^ fields are negligible for 
all narrow-band detected objects, but are as large as a factor of, say, ~ 2 — 5, for LAEs. 
These large differences probably come from field-to-field variation as well as Poisson errors. 
We evaluate the field-to-field variation in our survey area, cr^, with 



~ (^g-.lFovi'^DM/'^DM-.lFov) (4) 
^llFoV = [{{^,:lFoV-^,?)-^,]m (5) 

where a dm and ctdm-.ifoV are the rms fluctuation of dark matter in all the survey volume 
of 1 deg^ and Suprime-Cam's 1-FoV volume (~ 0.2 deg^), respectively. We calculate the 
flu ctuations of dark matter with the power spectrum, adopting the transfer function given 
by lBardeen. Bond. Kaiser. &: Szalayl (Il986l ). CFg^ipoV is the fluctuation of number density of 
LAEs for the 1 FoV. T,g,iFov and are LAE's surface densities of the 1 FoV and the average 
of survey area. 

Since these estimates of field-to-field variation are based on a large but single contiguous 
field, it is important to check whether our field is located at the sky of an overdense or 
underdense region. In fact, a large-scale overdensity or underdensity of Lya sources could 
also be produced by an inhomogeneous distribution of Lya absorbers (i.e. neutral hydrogen) 
along the line of sight. 

In Figure [11] we compare the surface densities of our LA Es with those on compl etely 
different s ky are as but selected with the same A^i?816 filter by lShimasaku et al.l (120061 ) and 



Hu et al.l (120041). We find tha t the surfa ce densities of o ur z = b.7 LAEs are consistent 



with those of IShimasaku et al.l (120061 ) and IHu et al.l (12004 ) within the scatters and Poisson 
errors of our 5 sub-fields. Moreover, Figure [TT] shows that Shimasaku et al.'s and Hu et al.'s 
measurements scatter around our average surface- density curve. Thus, we conclude that our 
single contiguous 1 deg^ field has no signature of overdensity or underdensity, and that our 
5 sub-fields represent the average field-to-field variation. 

In the same manner as the calculations of surface densities. Figure [12] plots color dis- 
tributions with the uncertainties of field variance for our LAE samples. In Section 15. H we 
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compare these color distributions (as well as the surface densities) with those reproduced by 
our Monte-Carlo simulations, and determine the best-estimates of luminosity function by 
fitting. 

3.2. Completeness and Contamination of the Samples 

First, we estimate the detection completeness of each narrow-band images, fdeti & 
function of narrow-band magnitude. We distribute 7875 artificial objects that mimic LAEs 
on our original 1 deg^ images after adding photon noise, and detect them in the same manner 
as for the detection of our LAE catalogs with SExtractor. We repeat this process 10 times, 
and compute the detection fraction. We find that the detection completeness is typically 
> 80 — 90% for relatively luminous sources which are 0.75 magnitude brighter than the 
magnitude cuts of our samples (i.e. A^S503 < 24.55, NB570 < 23.95, and NB816 < 25.25). 
The detection completeness is > 50 — 60% even for the faintest magnitude bin of our samples 
(i.e. NB503 = 24.8 - 25.3, iV5570 = 24.2 - 24.7, and NB816 = 25.5 - 26.0). 

Second, we estimate the contamination of our photometric LAE samples, which we need 
to consider in calculating the LP. We use our spectroscopic catalog of 3,233 objects (Section 
12. 4p . We define the contamination fraction, /cont, with 

/,o„t = 1-NI\JN^, (6) 

where A^lae ^au numbers of spectroscopically-identified LAEs and all spectro- 

scopic objects, respectively, lying in our color criteria (i.e. eqs[I]|3D. 

Since we have spectroscopic objects with no identification (see Section [X4|) . we calculate 
/cont for the following two extreme cases. If we omit the unclear objects, we find A^LAE/-^an — 
(29/29,14/14,11/11) for z = (3.1,3.7,5.7). Note that we see no obvious interlopers in 
our photometric samples. In case that all the no-identification objects are interlopers, 
^lae/^S = (7/8,6/7,3/4), where we take magnitude cuts of {NB5Q3, NB570, NB816) = 
(24.0, 24.0, 24.5) that are bright enough to be completely identified by our spectroscopy. 
Thus, the contamination rates are taken within the ranges of /^nt = (0 — 13%, — 14%, — 
25%) for z = (3.1,3.7,5.7) photometric LAE samples, respectively. Since these contamina- 
tion rates are negligibly small, we do not correct them for contamination. It is notable that 
statistical errors from contamination cannot be as large as ~ 30%. 



4. AGN and Multi- Wavelength Detection in our Samples 



We identify AGNs in LAEs witli our spectra and SXDS multi-wavelengtli images of 
radio, sub-mm, mid- infrared, and X-ray, and investigate AGN activities in out LAEs. We 
summarize the properties of these AGNs in Table HI In this section, we formally define the 
positional uncertainty of optically-identified LAEs, cXopt = FWHM/2.?,b, with a FWHM of 
narrow-band morphology. 



4.1. Spectra with AGN signatures 

We find 3 and 1 objects at 2; = 3.1 and 3.7 whose spectrum shows high-ionization lines, 
such as Nv, Civ, Hell, and Cm], originated from AGN activities. Figure [T^ plots spectra 
of these four objects. All of these AGNs show broad lines in these high-ionization lines 
with a line width of fpwHM ^ 1000km although the line width for one of these object, 
NB503-N-35820, has a large uncertainty. Radio bright objects, i.e. radio galaxies, are also 
included in our sample. The properties of these AGNs are presented in Table HI 

We also find that NB503-N-42377 shows a marginal (~ 2a) CiV emission. Although 



a typical high-z AGN has a line ratio of Lya/Civ= 5-10 (e.g. lMcCarthvlll993h . NB503- 
N-42377 has Lya/Civ> 25. With the facts of the marginal Civ detection and the large 
flux ratio of Lya/Civ, we rule out an AGN being the dominant power source of Lya for 
NB503-N-42377, and do not classify this object as an AGN. 



4.2. X-ray Detection 



We use a deep XMM-Newton EPIC image of SXDS, and identify an X-ray counterpart 
of LAEs with the combined X-ray catalog of SXDS ver 5.0 (Ueda et al. in prep.). This X-ray 
catalog is made of sources with a detection likelihood greater than 7 in either of 6 bands, 
0.3-0.5, 0.5-2, 2-4.5, 4.5-10, 0.5-4.5, and 2-10 keV (see Ueda et al. in prep). We calculate a 
combined positional error by taking a sum of square errors of optical and X-ray positions. 
Although there are 3 X-ray counterpart candidates with a < 2(T level of the combined error, 
we find by our visual inspection that one of the candidate, NB503-N-87126, clearly appears 
to be confused by the neighboring foreground objects. This X-ray counterpart is regarded as 
a foreground object. We thus identify 2 X-ray counterparts; 1 and 1 counterparts at 2; = 3.1 
and 3.7, respectively (Table H]). No counterparts are found for z = 5.7 LAEs. 
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4.3. Infrared Detection 



The Spitzer MIPS 24/xm-band image is investigated for infrared counterparts of LAEs. 
We use relatively sh allow data of the Spitzer Wide-area InfraRed Extragalactic (SWIRE; 
Lonsdale et al.l |2003| ) survey. The SWIRE MIPS data cover almost entire SXDS field of 
a 0.9 deg^, and miss only 13, 6, and 15 LAEs at 2; = 3.1, 3.7, and 5.7, respectively, 
which are located at the edge of our SXDS-W field (see Figure 1 of Morokuma et al. 2007 
submitted). The depth of MIPS data is cr = 48/iJy (Shupe et al. in preparation; see 
Ivison et al.l 120071 ). We run SExtractor to detect and measure source fiuxes with a 12"- 
diameter aperture. We apply the aperture correction of a factor of 1.698 to the aperture 
fluxes (see http://ssc.spitzer.caltech.edu/mips/apercorr/). We have cross matched 
our LAEs with MIPS sources detected at a > 5o" level, i.e. > 240/iJy. We flnd 3 (and 0) 
MIPS counterparts aX z = 3.1 (at z = 3.7 and 5.7) within the 3a positional allowance of our 
LAEs (Table i]). 



4.4. Sub-mm Detection 



We compare the positions of LAEs with SCUBA 850/xm sources of Submm HAlf-Degree 
Extragalactic Survey (SHADES: ICoppin et al.ll2006l ). The rms noise level is ~ 2.2mJy. The 
SHADES image covers an ~ 360 arcmin^ area at the center of SXDS-C, which corresponds 
to about 10% of our 1 deg^ SXDS field. In this SHADES region, we have about 16, 6, and 
3 LAEs ai z = 3.1, 3.7, and 5.7, respectively. These LAEs include one X-ray emitting LAE 
of NB503-C-49497 at z = 3.128 (see Tabl e HI) . The beam si ze of the sub-mm image is as 
large as FWHM = 14.7 arcsec. Following llvison et al.l (120071 ). we search counterparts with 
a positional error circle of radius 8 arcsec. We find no submm counterparts in our LAEs 
(including the X-ray emitting LAE) located at the SHADES area. We extend our search 
radius up to 12.5 arcsec for completeness. We find that a 2 = 5.7 LAE of NB816-C-90169 falls 



on the large search radius of a submm source of Sxdf850.5 (ICoppin et al.l 120061 ). However, 
on our optical images this submm source is very likely centered at a foreground galaxy of 
SXDS-iC-086385 with i' = 21.11. Thus, we conclude that none of our LAEs has a sub-mm 
emission down to ^ 3 — 5mJy in the SHADS area. 



^ All of X-ray and radio emitting LAEs listed in Tableware included in the MIPS area. 
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4.5. Radio Detection 



W e investigate radio pr operties of LAEs with deep VLA 1.4-GHz radio image of the 
SXDS (ISimpson et al.ll2006a] ). This radio image reaches an rms noise level of 12/xJy beam~^. 
We search a radio counterpart of LAEs in a 3x3 pixel box of the VLA image (1 pixel=1.25 
arcsec) . A radio source with a signal-to - noise ra tio of > 5 are identified as a radio counterpart 
in the same manner as ISimpson et al.l (j2006bl ) . We then visually inspect LAEs with radio 
contours, and find that two LAEs are detected in the radio image (Table H. 



4.6. LAEs Hosting AGN 



As summarized in Table HI we find that A z = 3.1 and 1 z = 3.7 LAEs have a detection 
in the multi- wavelength data. Figure [H] presents the spectral energy distributio n (SEP) of 



these AGN-LAEs, together with those of the av erage radio quiet /loud quasars (jElvis et al. 
1994h . We also plot a UV-optical spectru m of iTelfer et al.l J2002h : iRichards et al.l J2003h 



19941 ). We refer to these average 



that is normalized to the average SED of (jElvis et al. 
SEDs+spectrum as QSO templates. 

Since the detection limits of infrared and X-ray (and radio) bands are as large as 
uLiy > 10^^ (10^^) erg s~^, LAEs with any multi-wavelength counterparts likely host AGN. In 
fact, we have 3 spectroscopically-identified LAEs with multi- wavelength counterparts (Table 
H]). The spectroscopic classifications of these LAEs indicate that all of these LAEs have 
AGN. Thus, we classify all of 5(= 4 + 1) LAEs with multi-wavelength detection as AGN. In 
addition to these 5 AGNs identified by multi-wavelength data, we have 4 spectroscopically- 
identified AGNs. Since 3 AGNs have both a multi-wavelength detection and spectroscopic 
identification, the number of our LAEs with an AGN (hereafter AGN-LAEs) is 6 in total: 
4 and 2 at z = 3.1 and 3.7, respectively. No AGN-LAEs are found aX z = 5.7. Note that 
3 2; = 3.1 and 1 z = 3.7 AGN-LAEs out of 6 AGN-LAEs satisfy our color criteria of pho- 
tometric samples (eq. [1] or [21 see Table H]). These numbers (i.e. 3 and 1) are sufficiently 
smaller than those of entire photometric samples (356 and 101 for z = 3.1 and 3.7). Since 
the AGN fraction in our photometric sample is negligibly small, [f| AGN-LAEs do not largely 
contribute to the statistical properties of LAEs in the following sections. 



^ We find that a radio source is located around NB503-S-29853. Since the flux center of this radio source 
is obviously placed outside of NB503-S-29853, this radio source is likely associated with the other source(s). 

We discuss the AGN fraction of our photometric samples in Section 17.21 We address star-forming 
activities of LAEs with the multi- wavelength data in Section 17.31 



- 15 - 



5. Luminosity Functions 

5.1. Lya Luminosity Functions 

We derive luminosity functions (LPs) of LAEs a.t z = 3.1, 3.7, and 5.7 from our 
photometric samples. First, we calcul ate LFs with a s i mple classical naethod that was 
taken by most of previous studies (e.g. 



Ouchi et al. 


2003: 


Aiiki et al. 


2003; 


Hu et al. 


2004 



Malhotra fc Rhoadsl I2004J . We obtain the number densities of LAEs in each magnitude bin 
by simply dividing the observed number counts of LAEs in a given narrow band by the effec- 
tive survey volume defined as the FWHM of the bandpass times the area of the survey. Here, 
we calculate the Lya luminosity of each object with the response curves of narrow and broad 
bands by subtracting the continuum emission measured from the continuum magnitude from 
the total luminosity in the narrow band, assuming that Lya enters in the central wavelength 
of the narrow band. In this calculation, we use total magnitude of narrow-band images. The 
continuum emission is estimated by the narrow-band excess color (i.e. color of narrow-band 
and broad-band) defined with a 2" aperture, so as to keep high signal-to-noise ratios and 
to avoid the effects of source confusion on broad-band images with high source density. To 
check the accuracy of this calculation. Figure [15] compares Lya fluxes measured from our 
images and spectra for our spectroscopic samples. Figure [T5] shows that both measurements 
agree well within error bars for most of LAEs. These values mostly scatter within the range 
of difference of a factor of 2. 

Figures [TMTS] present the LFs from this classical method. Note that we corrected the 
detection-completeness by weighting with f^et measured in Section [221 To check the field-to- 
field variation and the accuracy of our results, we plot the estimates of LFs from the entire 1 
deg^ field (filled circles), together with those from the five ~ 0.2 deg^ fields (open symbols). 
In Figures [TMTTt we calculate the errors of field-to-field variations with eq. (jl]), and include 
these errors in the error bars of LFs for the entire 1 deg^ field. We find that the field-to-field 
variations are as large as a factor of ~ 2 — 5 among the five 0.2 deg ^ fields, although the 
typical scatters of the 0.2 deg ^ results are not far beyond the errors of Poisson statistics. 



Then we obtain the best-fit Schechter function (jSchechteii Il976l ) defined by 



(PiL)dL = (t)*{L/L*Y exp{-L/L*)d{L/L*). (7) 

The classical method is accurate when the narrow-band filter has an ideal boxcar shape. 
However, since the shapes of actual filters used in LAE surveys are rather close to a trian- 
gle, the classical method potentially suffers from the following two uncertainties, (a) The 
narrow-band magnitude of LAEs of a fixed Lya luminosity varies largely as a function of 
redshift. Thus Lya luminosity may be over- or under-estimated for some LAEs. (b) The se- 
lection function of LAEs in terms of EW also changes with redshift; the minimum EW value 
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corresponding to a given (fixed) narrow-band excess, such as i' — A^i?816 for z = 5.7 LAEs, 
becomes larger when the redshift of the object goes away from the redshift corresponding to 
the center of bandpass. 

In order to avoid such uncertainties, we perform Monte Carlo simulations to find the 
be st-fit Schechter pa r amet ers for the Lya LF. These Monte Carlo simulations are developed 



by IShimasaku et al.l (120061 ). and the details of simulations are described in their paper. We 
briefly summarize the procedure of these simulations. We generate a mock catalog of LAEs 
with a set of the Schechter parameters {a, (p*, L*) and a Gaussian sigma of probably dis- 
tribution of EW, (Tew- We uniformly distribute them in comoving space over the redshift 
range of all LAEs that can be selected with our narrow-band filter. We 'observe' these LAEs 
with the narrow- and broad-bands same to the real band responses (Figured]), and add to 
their flux densities photon noise corresponding to the actual observation. We select LAEs 
by the same criteria as for selecting the actual LAEs, and derive the number densities and 
color distributions from the mock catalog. We compare these results with the observational 
measurements of number densities and color distributions. By performing this set of simu- 
lations over a wide range of the parameters, we find the best-fit Schechter parameters by 
fitting. We compare the results of these simulations with the observational measurements 
corrected for the detection completeness, since the simulations are not affected by the de- 
tection completeness, fdet- In Figures [9lfT2| we plot the predicted surface densities and color 
distributions of LAEs for some parameter sets. 

The difference in for values is found to be insignificant. We fix the a value to -1.0, 
-1.5, and -2.0, and carry out fitting, varying 0*, L*, and ctew- Figures iQlfTT] present that 
all three models of a = —1.0, —1.5, and —2.0 reproduce the observed counts well, although 
shallower a gives a slightly better fit. This suggests that the Lya LF of LAEs is approximated 
well by the Schechter function. The best-fit parameters, thus obtained, are summarized 
in Table O W e adopt the a = —1 .5 results as the fiducial set of the best-fit Schechter 



parameters (see iGronwall et al.ll2007l ). Our best-fit Schechter parameters with a = —1.5 are 
0* = {9.2tll3Atll7.7tli) X 10-^ Mpc'^ and Ll^^ = (5.8^°;?, 10.2lii 6.8l|?) x 10^^ ^^g 
s-i at z = (3.1,3.7,5.7). 



The best-fit Schechter functions for a = —1.5 obtained above are shown in Figures [ 
with the solid lines. Note that the LFs from the classical method (filled circles) are consistent 
with the best-fit Schechter functions from our simulations. This means that the simple 
method gives a good approximation. This is probably because the uncertainties (a) and (b) 
in the classical method are negligible and/or cancel with each other. The number densities 
of the best-fit Schechter functions are slightly larger than those of the classical method. This 
slight difference can be explained by the underestimates of the number density and Lya 
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luminosity with the classical method due to the large survey volume of the top-hat redshift 
distributions and the small Lya flux measured with a bandpass of the triangle-shaped narrow 
bands. 

Note that these best-fit LFs are based on our LAE samples selected with a limit of EW 
(Table [3]). If we extrapolate the EW distribution down to EW = oA with our simulations, 
we can obtain estimates of the LFs for all LAEs with EW > OA. We refer a 0* value of 
LFs for all {EW > OA) LAEs as for 0o, and present the best-fit 0o Table El These 0o 
values only differ from 0* by ~ 10%. There are two reasons for this small difference: (i) 
The typical distribution of EW, o"eW) is much wider than that of EW limit of our selection, 
EWiim (see Tables |3] and [5]) . (ii) the EW range missed by our selection (0 < EW < EWum) 
is significantly shorter than that of our selection {EW > EWum) ■ The combination of these 
facts provide this small (~ 10%-level) difference between 0* and (pQ. 



5.1.1. Comparison with Previous Measurements 



We overplot previous measurements of Lya LFs for 2; ~ 3, 4 and 6 in Figure [191 Our 
LFs of 2; ~ 3 and 4 agree very well with all of the previous measure ments typically within 
the error bars. It is noticeable that the best-fit Schechter function of iGronwall et al.l (120071 ) 
(from the 0.3 deg^ survey) shows an extremely good agreement with ours within a factor of 
~ 0.1, although the luminosity range of their measurements is limited to the relatively faint 
luminosity (i.e. logL ~ 42 — 43 erg s~^). 

In the bottom panel of Figure (THl the previous measurements of 2; ~ 6 LAEs show 
lar ge scatters over P oisson errors. However, our measurements are consistent with those 



of IShimasaku et al.l (120061 ) (red open circles) as well as the recalculated LF of iHu et al. 
( I2OO4J ) (squares), both of which are based on deep and moderately- wide (~ 0.2 deg^) field 
imaging and spe c trosc opy. Moreover, z = 5.7 LF of the recent 2-deg^ imaging survey of 



Murayama et al.l (120071 ) (red open diamo nds) agrees with ours wit hin la error bars at log L > 
43.0 erg s~^ where t he measurements o f iMuravama et al.l (120071 ) are complete. (Note that 



and that the data points of iMurayama et al. 



the imaging data of Muravama et al. (|2007[ 



The best-fit Schechter function of 



lower than our LFs. Shimasaku et al 



are sh allower than ours by ~ 1 magnitude, 
J2007h are not corrected for incompleteness.) 



Ma 



20^ and 



hotra fc Rhoadd (120041) (dotte d line) is significantly 



Tapken et al.l (120061 ) have also pointed out 
that the measurements of iMalhotra fc Rhoadd (|2004f) are lower than their L Fs. Our z = 
5.7 LF is consistent with the measurements of lAjiki et al.l (120061 ) (stars) and iTapken et al. 



( 120061 ) (hexagon) within the error bars, but our measurements would be systematically higher 
than theirs. This difference is probably originated from i) the method of Lya luminosity 
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estimates from photometric data, ii) the detect ion-complet e ness c orrec tion and iiii) f i eld-to - 
field variation raised by smaller field surveys of lAjiki et al.l (120061 ) and lTapken et al.l (120061 ) . 
The differences in selection criteria for LAEs between our sample and the others are probably 
not a main reason for this discrepancy, since, for our data, changing the i' — NB816 threshold 
over 1.2 - 1.5 only makes a difference in number density by 11%. 



5.1.2. Evolution of Lya Luminosity Function 

Our LFs are appropriate to identify an evolutionary trend of LF, since our LPs from 
z = 3.1 to 5.7 are derived by the same procedure with similar data sets from the same 
instrument. In Figure [T9| we compare our Lya LFs of LAEs at 2; ~ 3, 4, and 6. The LFs do 
not change within error bars at these redshifts. The left panel of Figure [20] shows the error 
ellipses of Schechter parameters of our LFs at 2; = 3.1 — 5.7. All the 2(T-error contours well 
overlap each other. At the fixed slope of a = —1.5, the best-fit parameters of LFs change 
from z = 3.1 — 5.7 only by factors of 1.8 and 2.7 in L* and 0*, respectively. Since these error 
ellipses are obtained by our observational data whose limiting Lya EW is different at each 
redshift, we also plot in the right panel of Figure [201 the ellipses for all LAEs with a positive 
emission {EW > 0) which are obtained from the extrapolation with our simulations. In this 
case, the overlaps of z = 3.1 and z = 3.7 error ellipses become smaller. However, the error 
ellipses of z = 5.7 LAEs are much larger than those of z = 3.1 and 3.7. These large z = 5.7 
error ellipses still appear to be consistent with z = 3.1 and 3.7 error ellipses. The results 
for LAEs with EW > are similar to those for the observed LAEs. This is because the 
number density of observed LAEs (0*) differs from that of all LAEs (0o) only by ~ 10%, as 
discussed in Section 15. 1[ Thus, we conclude that Lya LFs of LAEs do not evolve by more 
than a factor of 2-3 either in luminosity or number density, and that drastic evolution in the 
Lya LF of LAEs from z ~ 3.1 to z = 5.7 is ruled out. It should be noted that the central 
positions of z = 3.7 contours appear to be shifted to brighter L* and lower 0* than those of 
z = 3.1 and 5.7, although this shift is not significant within 1 — 2a levels This small shift 
may indicate a moderate evolution of LFs between these redshifts. However, it is possible 
that AGNs at the bright-end LF more strongly affect the Schechter fit of z = 3.7 LF than 
that of z = 3.1 and 5.7 LFs, due to the shallow detection limit of our z = 3.7 LAEs (see 
Section [T6l) . 

Note that Lya fiuxes from high- 2; objects are generally attenuated by neutral hydrogen 
of intervening inter-galactic medium (IGM). It is known that high- 2; galaxies have an asym- 
metric Lya emission whose blue si de of line is more strongly absorbed by neutral hydrogen 



of IGM than low-z galaxies, (e.g. [Hu et al.|[2004j : [Shimasaku et al.[ [20061 ). Thus, the Lya 
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la 



LFs in Figure [19] (solid lines) should be called apparent-Lya LPs. Due to the absorption of 
IGM, we cannot directly measure intrinsic-Lja fluxes that are emitted from LAEs. How- 
ever, we can constrain the evolution of intrinsic-Lja LFs with apparent-Lja LFs. First, the 
lower limits of intrinsic-Lja LFs are obviously given by the apparent-Lya LFs (solid lines of 
Figure [T9|) . Second, we can estimate the intrinsic-Lya luminosity with the average Gunn- 
Peterson (CP) optical depths. The GP optical depths a r e esti mated in variou s studies . Here , 
based on the results of IPan et al.l J2006h fPan). iMadaul Jl995h (Madau), andlMeiksnJ fl2006h 
(Meiksin), we calculate the ratios of intrinsic Lya flux, f^^^, to apparent Lya flux, f^y^. 
Assuming that IGM absorbs a blue half of symmetric Lya emission line, we estimate the 
ratios fl^^/ flyi with the formulae of (Fan, Madau, Meiskin) to be (0.85, 0.81, 0.84), (0.76 
0.73, 0.76), and (0.53,0.54,0.52) at z = 3.1, 3.7, and 5.7, respectively. Since the f^/f/" 
ratios change only by 3 — 5% between these three estimates, we use the classic lMadaulll995l 's 
model. The dashed lines in Figure UTl correspond to the intrinsic-Lya LFs that we estimate. 
This figure indicates that intrinsic-Lya LFs may evolve significantly, and that luminosities 
and/or number densities of LAEs are intrinsically brighter/higher at 2; = 5.7 than z = 3.1. 

Note that we estimate these intrinsic-Lya LFs with two assumptions. The first as- 
sumption is that the density of neutral hydrogen atoms in the intervening IGM is the cosmic 
average value at each epoch. However, hydrogen clouds around LAEs are likely more ionized 
by the UV radiation from LAEs. On the other hand, the chance of intervening galaxies 
with molecular hydrogen probably increases around LAEs, due to clustering of galaxies. 
The second assumption is that the intrinsic-Lya emis sion of LAEs has a sy mmetric profile 
whose central wavelength is not red- or blue-shifted. IShapley et al.l (120031 ) find in a com- 
posite spectrum of 2; = 3 LBGs that Lya emission is redshifted by 360 km s~^ from the 
average of continuum absorption lines. Because the typical Lya FWHM of LAEs is as small 
as < 300 km s"^ (e.g. iRhoads et al1l2003l : IHu et al.lliooi : Ishimasaku et al.ll2006h . the IGM 
will absorb only a small fraction of Lya photons if they are redshifted by about 360 km s~^. 
However, this is true only if Lya properties of LAEs are similar to those of 2; = 3 LBGs. In 
fact , the typical Lya F WHM of LAEs is smaller than that oi z = 3 LBGs (450 ± 150 km 



s~^ IShapley et al.l 120031 ). It is possible that the velocity offset of Lya emission may also be 
smaller in LAEs than in LBGs (360 km s~^). In this way, the inferred intrinsic-Lya features 
(such as those found in Figure [?I1) are real only if these two assumptions are correct. We need 
high-resolution spectroscopy for LAEs and sophisticated numerical simulations to assess the 
intrinsic-Lya LP, and to derive a reliable conclusion on its evolution. 

In summary, it is interesting that apparent LP of LAEs does not show a significant 
evolution over the long redshift range of z = 3.1 to 5.7. We point out the possibility that 
evolution of IGM and LAEs may be tuned to realize this no evolution aX z = 3.1 — 5.7, if 
the amount of external absorption changes with redshift. 



- 20 - 



5.2. UV Luminosity Functions 



We derive the UV LFs of our LAEs at z = 3.1, 3.7, and 5.7. We use our LAE catalogs, 
and obtain UV-continuum luminosities from magnitudes of a broad band whose bandpass 
does not include Lya but only continuum emission. We choose R magnitudes for z = 3.1 and 
3.7 LAEs, and z' magnitudes for z = 5.7 LAEs, where the corresponding wavelengths are 
~ 1300 — 1600A in rest frame. We apply no k-corrections to the UV magnitudes, since the 
difference of UV magnitudes in this narrow-wavelength range is negligibly small for LAEs 
with a fiat continuum (Section 16. 3p . We calculate number densities from the survey volume 
and the number counts as a function of UV magnitudes by weighting with the detection 
completeness (Section 13. 2p . Note that we need a careful estimation of total magnitudes of 
broad band. Since the aperture sizes and shapes of MAG_AUTO magnitudes are defined by 
the sources on the narrow-band images, the MAG_AUTO magnitudes of broad-band images 
could be contaminated by the neighboring continuum sources. Thus, we identify and deblend 
sources on our broad-band images with SExtractor, and cross-match the broad-band sources 
with our LAEs. If a LAE do not have a corresponding broad-band source detected over a 
3a level, we just use photometry of a 2"-diameter aperture. We refer to these magnitudes 
as a total magnitudes of broad band, and derive UV LFs with these magnitudes. For 
comparison, we also calculate the LFs from the simple 2"-diameter aperture photometry 
with no completeness correction. These simple estimates provide lower limits of UV LFs. 



Figure [22] shows the UV-LFs of our LAEs, together with the previous measurements of 
L AEs and dropout galax ies (or LBGs) at e a ch red shift. Our UV-LFs are consistent with those 
of IHu et al.l (120041 ) and IShimasaku et al.l (120061 ) at z = 5.7, although their measurements 
show slightly higher number densiti es than ours. To check the consistency, we calculate 
z = 5.7 LF in the same manner as IShimasaku et al.l (120061 ) who simply used 2"-diameter 
aperture magnitudes with aperture co r rectio ns estimated with their NB816 image. Our 
LF by the method of IShimasaku et al.l (120061 ) presents a very good agreement with that of 
Shimasaku et al. ( 20061 ). However, we find that the aperture corrections defined with NB816 
band overestimate the broad-band magnitudes. This is because the Lya morphology in the 
narrow band is generally more extended than UV morphology on the broad band. Thus, we 
take our best estimates based on the MAG_AUTO magnitudes with the source deblending 
on the broad-band images. 



In Figure [221 we find that UV-LFs of LAEs do not change at 2; ~ 3 — 4, but evolve from 
2; ~ 4 to 5.7. There exist an emergence of UV-bright LAEs at z = 5.7. We compare these UV 
LFs with those of dropout galaxies. At z = 3—4, number densities of our LAEs are as much as 



10% of dropout gal a xies down to M\ 



20, which corresponds to ^ M*jy or ^ M 



(ISteidel et al.l Il999l : lOuchi et al.l l2004al : ISawicki fc Thompson! l2006l : iBeckwith et al. 



+ 1 



2006 
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Yoshida et alJbood : IVanzella et al.ll2006bf ). Note that the EW hmit of our z = 3 - 4 LAE 



samples is ~ 40 — 60A. This ratio of Lyo; em itting galaxies to dropout galaxies is consistent 
with that reported by IShapley et al.l (120031 ) who find that 25% and ~ 2% of their z = 3 
dropout galaxies have Lya emission line with EWq > 20A and EWq > lOOA, respectively. 
On the other hand, at z = 5.7, the number densities of our LAEs seems comparable to, or 
at least more than 50% of, dropout galaxies at z ^ 6, where the EW li mit of our z = 5.7 
LAE is about ~ 30A. Spectroscopic studies of iDow-Hygelund et al.l (120061 ) report that 6 Lya 
emitting galaxies with EWq > 20A are found in their i'-dropout sample of CL1252 and HUDF 
parallel fields, and conclude that about ~ 30% of z'-dropout galaxies have a Lyg emission 
line. T his result is consistent with the meas urements of preyious studies (IStanway et al. 



2004al lbl: IVanzella et"alll2006ah . More recently, Istanwav etaP ( 120071 ) show that 16 out of 24 
i'-dropout galaxies do not have a Lya emission line with EWo > 25 A in the GLARE/HUDF 
field. Most of the previous spectroscopic studies of z'-dropout galaxies suggest that ~ 30% 
of i'-dropout galaxies down to ~ L* have a Lya emission line with EWq > 20A. It is not 
clear why our ratio of dropout galaxies and LAEs is higher than those of spectroscopic 
results. There are three possible explanations. The first is the difference of the definition 
in t he ratio. We compare the rati o of LAEs to dropou t galaxies at the same redshift, 
but iDow-Hygelund et al.l ( l2006l ) and IStanway et al.l (120071 ) obtain the ratio of Lya emitting 
dropouts to all the dropouts. If all the LAEs are identified by dropout technique, these 
two ratios should be comparable. However, LAEs with no significant continuum are likely 
missed by the dropout selection. Moreover, the i dropout selection systematically misses 
the strong Lya emission lines at z < 6. The selection differences may cause the discrepancy 
between these two ratios e stimated by the different methods. In fact, the theoretical model of 



Dijkstra fc Wyithd (120071 ) suggests that LAE selections identify more population III galaxies 



with a large EW than dropout selections, and that the EW distribution of LAEs should not 
be the same as that of dropouts. The second explanation is that LPs o f dropout galaxies have 
systematic errors. Indeed, we compare the LPs of dropout galaxies of iBouwens et al.l (120061 ) 
to estimate the ratio. However, the measurements of dropouts' LPs scatter by a factor of 2 
or more in the different estimates by various authors. We show the uncertainties of estimates 
with the gray shade in Pigu re [221 The third explanatio n is the evolution of UV LPs from 
z = 5.5 to 6.5 as claimed by iDow-Hygelund et al.l (120061 ). The UV -LPs of dropout galaxies 



decrease in number density and/or lumino s ity from z^Ato 6 (e.g. iLehnert &: Bremerll2003 



Ouchi et al. 2004a; Shimasaku et al. 2005; Bunker et al. 



2004J'). and even more significantly 



evolve from z ~ 6 to higher redshifts of ^ ^ 7 — 8 (jBouwens et al. 2006). Since i-dropout 
galaxy samples distribute over z = 5.5 — 6.5, the LPs of i-dropout galaxies are measured at 
higher redshifts in average than that of our LAEs (i.e. z = 5.7). This would not allow us to 
obtain the accurate ratio of Lja emitting population to dropout galaxies. Por these reasons, 
it is likely that the simple comparison between LPs of dropout galaxies and LAEs would not 
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give a ratio between all z ~ 6 objects and Lya emitting population. 

It is interesting to compare the shapes of UV LF of LAEs and dropout galaxies. At 
2; = 3 — 4, we find a significant deficit of UV-brig ht LAEs with respe ct to dropout galaxies. 
This deficit of UV-bright LAEs are reported by lOuchi et al.l (120031 ) at z = 4.9 based on 
the comparison of dropout and LAE LPs. This te n dency is also cons i stent with the results 
of Secti on 16.21 and the claims of lAndo et al.l (120061 ) ; IShimasaku et al.l (120061 ) ; IVanzella et al. 
(l2006bl ) who find that UV-bright galaxies do not have a large Lya EW from spectroscopic 
samples at 2; = 5 — 6. On the other hand, this tendency cannot be clearly seen for z ~ 6 in 
Figure [221 probably due to the large uncertainties of LF of dropout galaxies. 



We fit the Schechter function, 



ij{M)dM = C<p* exp {-C{a + 1)(M - M*) - exp[-C(M - M*)]} c/M, 



(8) 



to the UV LFs, where C = 0.41n(10): a is the power-law slope, (j)* is the normalization 
factor, and M* is the characteristic absolute magnitude. We fit the Schechter function with 
the measurements over a 5a level (filled circles in Figure [22]) . Since the faintest points of our 
z = 3.7 LAEs appear to be suffered from the incompleteness of source detection, we do not 
use the point for the fitting. Since our LAE samples are not enough to derive a, we fix the 
slope to a = —1.5. For our z = 3.7 sample, we fix the M* to the best-fit M* value of our 
z = 3.1 LAEs. We summarize the results of the fitting in Table El 

We find no evolution at z ~ 3 — 4, even considering the error bars. The difference 
between z = 3.1 and 3.7 LFs is within 50% in number density or luminosity. On the other 
hand, UV LFs of LAEs evolve significantly from 2; ~ 3 — 4 to 6. If it comes from the pure 
luminosity evolution, the UV LFs brighten by a ~ 1 magnitude from 2 ~ 3 — 4 to z = 5.7. 
Note that the EW limits of our selections are different between these redshifts. However, 
the difference of sample selection only changes by ~ 10% in number density, according the 
results of our simulations (Section 15.11) . Thus, this evolutionary trend of UV LFs is not 
affected by the selection effects. 



5.3. Ly«- and UV-Luminosity Densities Contributed from Lya Emitters 

We calculate Lya and UV-luminosity densities of our LAEs from the luminosity func- 
tions (LFs) derived in sections 15.11 and 15.21 We show Lya and UV-luminosity densities in 
Tables and El respectively. These tables present Lya (UV) luminosity densities down to 
our observational detection limits, pLya (Puv)' total luminosity densities down to zero 
luminosity by extrapolation, pLya (Puv)- These observational and total luminosity densities 
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provide approximated lower and upper limits. For a reference, we estimate total Lya lumi- 



nosity densities of LAEs with EW > 0, from our simulations, which are shown under Pql* 
in Table E 



ya 



In Tabled the total Lya luminosity density is ~ 10^° erg^^s^^Mpc^^ at z = 3.1 — 5.7, 
and does not change within a la level from 2; = 3.1 to 5.7. Table [6] shows that the total 
UV luminosity density would increase from ~ 4 x 10^^ erg^^s~^Hz~^Mpc^^ to ~ 8 x 10^^ 
erg~^s~^Hz~^Mpc~^ at ^ = 3 — 4 to 5.7, although the difference is only a 1 a level because 
of large errors. This possible increase of UV luminosity density could be due to the fact that 
UV bright LAEs increase at 2; = 5.7 as shown in Figure 



Note that the entire Lya (and UV) emission is not only originated from star-formation 
activities, but also from cooling radiation, shock winds, and AGN, and that Lya (and UV) 
emission is attenuated by dust and IGM. However, we roughly estimate the contribution of 
star- formation from LAEs with those luminosity densities to check the consistency with the 
previous studies. We use the simple prescription of Lya luminosity, Li^ya, and star- formation 
rate, SFR, 

SFRiU^yr-') = U,^{eig s-i)/(l.l x 10^^)^ (9) 



combining the relation of Ha luminosity and star- formation rate (iKennicuttl Il998l ) and the 



case B appr oximation (Brocklehurstlll971l). With the assumptions of Salpeter IMF and solar 



metallicity (iMadau. Pozzetti. fc DickinsonI Il998l ). the relation between UV luminosity and 



star- formation rate, is: 

SFR(M0yr-i) = Luvi^^g s-^Hz-i)/(8 x 10^^), (10) 

where Luv is the UV luminosity measured in 1500A. In Figures [T^ and [221 we plot ticks of 
star- formation rates corresponding to Lya and UV luminosities based on these formulae of 
eqs ([9]) and ffTOl) . Assuming these equations, SFRs of a typical (L*) LAE are SFRilsAQji'^) = 
{5.3toi,9.3tii6.2tl-l) and (4.5 ± 2.0, ~ 4.5,9.4 ±7.0) at z = (3.1,3.7,5.7), which are 
estimated from Lya and UV-continuum L*, respectively. Thus, a typical LAE has a SFR of 
~ 5 — 10 Mq yr~^ at z = 3.1 — 5.7, although these SFRs are correct under the assumptions 
given by eqs. and ffTUl) . 

We estimate the nominal lower limits of cosmic SFR density (SFRD) from the luminosity 
densities down to our detection limits, PLy«- From Table [5], the observed Lya luminosity 
density, pg^„ = (4.8^1:^,2.41°;^, 3.61?;^) xlO^^ erg-^s-^Mpc'^ for redshifts of (3.1, 3.7, 5.7). 
If we naively estimate SFR densities (SFRD) from these Lya luminosity densities with eq. 

we obtain SFRD = (4.3lJ:^, '^■'^-ot 3.2+?;8) x lO'^ Moyr-^Mpc-^ for z = (3.1, 3.7, 5.7). 
Similarly, we estimate SFRDs from the observed UV luminosity densities, Puv? which are 
(1.2+J:^, 0.5+1]:}, 0.9lg:^)xl025 erg-^s-^Hz^iMpc-^ for redshifts of (3.1, 3.7, 5.7). Based 
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on eq. ffTOj) . we obtain SFRD 



1.5^8 ^,0.71°-!, 1.1 



0.1' 



+3.U 

-o.sJ 



X 10-3 Meyr-^Mpc-^ for z 



(3.1, 3.7, 5.7 ). These Lyg and UV SFRDs are consistent with those obtained by the previous 



studies (see iTaniguchi et aLll2005l and references therein). 



If we estimate the nominal total SFRDs from the total luminosity densities of p{^- 



,tot 



and that are given by the extrapolation of luminosity functions, we obtain SFRD = 
(8.5l}j, 5.71};^, 8.3+^:^) x lO'^ and (5.5l|^, 5.1+°;^, 9.41^6^) x IQ-^ Mgyr-^Mpc-^ for z = 
(3.1,3.7,5.7), respectively. We compare th e SFRD estima t ed fro m Lya LF of 2; = 3.1 
LAEs that is recently obtained by MUSYC. iGronwall et al.l (120071 ) have found that SFRD 
contributed hj z = 3.1 LAEs is ^ 6.5-8.6 x lO'^ Mgyr-^Mpc'^. Our result of 8.5;^! x lO'^ 
M0yr~^Mpc-3 is very consistent with their finding. 

The SFRD from UV luminosity density may increase from 2; = 3.1 to 5.7 by a factor of 2, 
but again the large errors do not distinguish between evolutionary effects and uncertainties. 
If we take into account the large error bars, SFRD ~ 5 — 9 x 10"^ MQyr-^Mpc-^ at 
2 = 3.1 — 5.7. We compare these SFRDs with those estimated from dropout galaxies. The 
results of SFRD measurements with no dust corre ction show that th e SFRDs of dropouts 



are 



2 - 3 X 10-2 Mq yr-i Mpc-^ at ;z ~ 3 



Hopkinsll2006l ). Thus, LAEs contribute roughly about 
at z = 3 — 6. 



6 fISteidel et al.lll999l : iBouwens et al.l 12006 
20 



-40% of the entire cosmic SFRD 



6. Properties of Lya Emitters 

6.1. Statistics of Lya Equivalent Width 

Figure [23] presents the histogram of rest-frame equivalent width (EWq) for our LAEs. 
We calculate EWq of our LAEs from our photometric measurements by the modeling de- 
scribed in Section 15.11 We then obtain two estimates of the EWq distribution by (i) taking 
the best measurements of EWo for each LAE and (ii) summing the probability distribution 
of EWo that is defined by errors of EWq. We show in Figure [23] these two estimates with 
(i) histograms (+Poisson errors) and (ii) solid curves. We find that the results of (i) and (ii) 
agree fairly well for z ~ 3 and 4 LAEs. However, for our z ~ 6 LAEs, we see a significant 
difference at EWo = 100 — 150A. This discrepancy between (i) and (ii) for z ~ 6 LAEs is 
probably due to the large errors of the relatively shallow off-band (i.e. 2;'-band) photometry. 
Thus, we regard the difference of (i) and (ii) as the uncertainties of our measurements. The 
gray histograms and curves indicate the EWq distribution for all of our photometrically- 
selected LAEs. Since our three LAE samples at z = 3.1, 3.7, 5.7 do not have the same limits 
of Lya luminosity and EWo, we make subsamples of LAEs that share the same detection 
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limits for comparison. We set logL(Lya) > 42.6 erg s"^ and EW'^^ > 70 - 80A for the 
limits of our subsamples, because our three LAE samples are complete in these limits. The 
histograms of these subsamples are plotted with black lines at each panel. We compare these 
histograms with those for spectroscopically identified LAEs. Although the statistical accu- 
racy is limited for the spectroscopic LAE samples, the EWq distributions of spectroscopic 
and photometric samples are consistent. 

In Figure [M] we plot gray areas showing the EWq distributions for the subsamples with 
the uncertainties of (i) and (ii). We find no clear evolution of EWq distribution from 2 = 3.1 
to 5.7 with the uncertainties of our measurements. Since these EWq distributions are based 
on the apparent EW, we correct for the absorption of IGM with the average attenuation at 
z = 3.1, 3.7, 5.7 (see Section 15.1.21) . We present in Figure [2H the distribution of intrinsic 
EWq (i.e. EWq corrected for IGM absorption) with cyan, blue, and red regions for z = 3.1, 
3.7, and 5.7, respectively. The IGM corrected histograms of ^ = 3.1 and 3.7 appear to be 
comparable within the large uncertainties of our measurements. The histogram of z = 5.7 
(red histogram) would be systematically different from those of ^ = 3.1 and 3.7, and could 
imply that z = 5.7 LAEs may tend to have a large intrinsic EWq. However, the uncertainty 
of EWq measurements is as large as the difference of histograms between at z = 3.7 and 
z = 5.7. We cannot reject the possibility of no evolution of intrinsic EWq from z = 3.7 to 
5.7. 

We investigate the fraction of large- EW LAEs with the subsamples. We define a thresh- 



old o i EWn'^^ > 240A with an intrinsic EW, following previous stud ies (e.g. iMalhotra &: Rhoads 



2OO2I : iDawson et al.ll2004l : IShimasaku et al.ll200a ISaito et al.ll2006f ) . This threshold of EW^ 



rint 



corresponds to the apparent EW of EWq^"^ = 192A, 175A, and 130A, assuming the average 
attenuation of IGM (Section I5.1.2p . We find that the fractions and errors of the large- EW 
LAEs in our subsamples are {z, fraction, fraction within 95th percentile)= (3.1, 0.21, 0.11- 
0.29), (3.7, 0.26, 0.16-0.40), and (5.7, 0.25, 0.15-0.33) by calculations of (i). From the 
calculations of (ii), we obtain the similar fractions of 0.20, 0.20, and 0.27 for z = 3.1, 3.7, 
and 5.7 LAEs. There is a possible but slight implication of the evolution of the fraction 
of large-EW LAEs from z = 3.7 to 5.7 in the results of (ii). However, again, this possible 
evolution is very uncertain due to the large errors of EW measurements. The best estimates 
of the fraction are 20 — 30% for z=3.1-5.7. Thus, we conclude with these uncertainties that 
the fraction of large-EW LAEs is 10-40% aX z = 3.1 — 5.7, which does not significantly change 
from 2; = 3.1 — 5.7 within this percent range. 



We compare our fraction of large-EW LAEs with those of IShimasaku et al.l (120061 ) for 



^ The ranges of 95th percentile are obtained by bootstrap samphng. 
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z = 5.7 LAEs. Since IShimasaku et al.l (120061 ) calculate this fraction with the threshold of 
apparent EWq of > 120A (c.f. 130A in above calculations), we use this threshold of apparent 
EWq > 120A here. With this threshold, we find the fraction and 95th-percentile range of 
0.34 and 0. 24-0.42 for our z = 5.7 LAEs. These numbers are vary consistent with those 
claimed by lShimasaku et al.l (120061 ) who obtain 30-40%. 



6.2. Relation between Lya and UV-continuum emission 



Figure[25lpresents rest-frame Lya EW (EWq) as a function of UV-continuum magnitude 
for our spectroscopic samples. The effective wavelength of UV magnitude is ~ 1300-1600A. 
The details of UV- magnitude measurements are presented in Section [5^ Figure [25] also plots 
EWq from previous studies of dropout galaxies and LAEs. On the right panel of Figure [251 
EWq of 2 ~ 6 objects show a clear defici t of large EWn ob jects in the bright magnitude 
(Muv < —2 1.5). This trend is fou n d bv lAndo et al. ( 2006 ) for z = 5 — 6 galaxies, and 
discussed in IShimasaku et al.l ( l2006l ): IStanway et al.l ( 120071 ) for z ~ 6 objects. We confirm 
their findings on our plot of 2; ~ 6 objects. For z ~ 3 and 4 objects shown in Figure [25], we 
find the trend very similar to that of 2; ~ 6. We see no large EWq objects wit h EWn > 80A 
for UV -bright (Myy < —21.5) objects. At z = 3, this trend is also reported by lShapley et al. 



(120031 ) who find that the Lya emission strength increases toward fainter magnitudes in their 
spectra of LBGs. 



6.3. UV-Continuum Color and Extinction 



We present UV-continuum colors as a function of UV magnitude for our spectroscopi- 
cally identified LAEs at 2; = 3.1 and 3.7 in Figure [261 This figure compares those of dropout 
galaxies at the similar redshift. Except for AGN, typical 2 = 3 — 4 LAEs have i' — z' (or 
R — i') colors of < 0.05. The UV-contin uum colors of our LAEs are consistent with those 
of 2 = 4.5 LAEs (IFinkelstein et al.l 120071 ). On the other hand, dropout galaxies (or LBGs) 
generally have colors of ~ 0.1. Thus, LAEs appear to be bluer than dropout galaxies in a 
UV-continuum color. This color difference indicates that LAEs are younger or less-dusty 
than dropout-galaxy population. We discuss more details in Section 17. 5[ 
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6.4. Composite Spectra 

We make a composite spectrum eacti for our z = 3.1, 3.7, and 5.7 LAEs. From our 
spectroscopic samples, we exclude AGNs (Table H]) and LAEs witli uncertain identification 
(Table [2]). We also do not include the spectrum of NB503-N-42377 that has a marginal Civ 
line, so as to obtain the average of typical LAE spectra (see Section H?T]) . Thus, we use 36, 25, 
and 15 spectra of z = 3.1, 3.7, and 5.7, respectively. Since some of spectra have systematic 
differences originated from a faint AGN and/or errors of cosmic rays and subtraction of sky 
etc., we calculate a mean flux with rejections of two largest and smallest values at each 
wavelength. Figure [27] shows the composite spectra of our LAEs whose average redshifts 
are (z) = 3.13, 3.68, and 5.68. We have also performed the median stacking, and find no 
qualitative differences between the rejected-mean and median results. The Lya luminosity 
of the composite spectrum is L{Lya) = (6.3 ± 0.1,8.5 ± 0.2,9.5 ± 0.3) x 10^^ erg s~^ at 
(z) = (3.13,3.68,5.68). We have a signal-to-noise ratio high enough to detect UV continua 
in (z) = 3.13 and 3.68 spectra, but not in (z) = 5.68 spectrum. 

The composite spectrum of (z) = 5.68 LAEs has Lya emission line with a clear asymme- 
try, which is a typical emission feature of high-2; objects. The spectra of (z) = 3.13 and 3.68 
do not show a significant asymmetry in Lja. This is because a typical spectral resolution of 
our z = 3.1 and 3.7 LAEs is not high enough to resolve the asymmetry of Lja {R ~ 500; c.f. 
R ~ 1000 for z = 5.7 LAEs). Moreover, the IGM absorption at ^ = 3 — 4 is not as strong as 
z = 5.7, which provides a weak asymmetric feature in Lya emission at 2; = 3 — 4. 

We find no significant absorption lines in these composite spectra, probably due to the 
combination of the weak metal absorption of LAEs and the relatively poor signal-to-noise 
ratios of our composite spectra. § The composite spectra of {z) = 3.13 and 3.68 cover 
the wavelength of HellA1640 emission. The inset panels of Figure [27] magnify the spectra 
around this wavelength. No significant feature is found for Hell at (z) = 3.13 and 3.68. 
We calculate upper limits of Hell flux, fueii, and estimate the 3a upper limits of the ratio, 
fueii/ fhya, whcrc /Lya is Lja flux. We obtain meaningful 3a upper limits of the ratio, 
fneii/ fhya = 0.02 and 0.06 for {z) = 3.13 and 3.68 LAEs, respectively Since Hell emission is 
an indicator of primordial populations, i.e. population HI stars and cooling radiation, our 
upper limits place constraints on primordial populations at z = 3 — 4. Al though theoretical 



models predict the existence of primordial populations at these redshifts (jScannapieco et al. 



20031 : 1 Jimenez fc HaimanI 120061 ). we have found no signature at the level of the upper limits 



of the present observations. Since the ratio of fueii/ fhya is predicted as srnall as ~ 0.1 — 0.00 1 



for population HI star formation and cooling radiation (e.g. lSchaererll2003l : lYang et al.ll2006l ). 



^ The Sill absorption of (z) = 3.13 LAE may exist, but the significance level is not high. 
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we need future surveys whose upper limits reach at the level of /neii/ZLya — 0.001. 



7. Discussion 

7.1. Evolution of Luminosity Functions 

We have found four pieces of evidence that provide a general picture in the evolution 
of LAEs: (1) apparent-Lja LPs of LAEs show no evolution from z = 3.1 to 5.7. (2) These 
apparent-Lja LPs and the average Gunn-Peterson optical depth indicate that intrinsic-Lya 
LPs of LAEs must have brighten in the luminosity and/or increased in the number density 
from z = 3.1 to 5.7. (3) The UV LPs of LAEs show the increase of luminosity and/or number 
densities from z = 3.1 to 5.7. (4) The rest-frame intrinsic-Lya EW, EW^'^^ , distribution of 
LAEs does not change or positively evolves from z = 3.1 to 5.7. 

The evolution of intrinsic-Lya LPs is caused by the emergence of Lya emitting objects 
with a bright UV continuum or the increasing trend of Lya EWo for LAEs at z ~ 5.7. There 
are three possible explanations for the evolution of LAEs. Pirst, the amount of dust mixed in 
HIT regions may decrease in average from 2; = 3.1 to 5.7, and the extinction of Lya luminosity 
through the resonance scattering decreases selectively. This implies that at 2; = 5.7 even UV 
bright galaxies have a detectable Lya emission line. The second explanation is that the young 
Universe at 2 = 5.7 has numerous exotic Lya populations such as metal-poor population III 
and/or cooling clouds. The third explanation is that the fraction of Lya-absorbed galaxies is 
lower at higher redshifts, due to an decrease in the amount of HI gas in and around galaxies 
with redshift. In fact, some z ^ 3 LAEs sh ow such an HI abs orption originating from 



an inflow and/or outflow in their spectra (e.g. iTapken et al.l 120071 ). although the evolution 



of such HI clouds has not been observationally understood. The CDM-based hi erarchical 



model predicts that galaxies at higher redshifts are generally more compact (e.g. iMo et al. 



19981 : ISomerville et al.ll2006l ). It is possible that the star formation in such compact high- 2; 
galaxies is so efficient that cold gas is immediately converted into stars. If high- 2; LAEs are 
such compact galaxies, they will have little HI gas and thus their Lya emission from HII 
regions can escape without absorption. If this hypothesis is true, we may be able to explain 
the evolutionary connection between LAEs and dropouts: the internal mass density of LAEs 
decreases with cosmic time by mergers and thus they evolve into dropout galaxies that are 
more massive and have more HI gas. 

It should be noted that the evolutionary trend of UV-LP of LAEs is opposite to that of 
dropout galaxies. The UV-LPs of dropout galaxies sh ow a significant decrease i n luminosity 



and/or number density from 2; ~ 4 to 2; ~ 6 (e.g. iLehnert &: Bremen l2003l : lOuchi et al. 
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2004al : IShimasaku et al.ll2005l : iBunker et al.l 12004 ). In case of lu minosity evolution, U V-LFs 
of dropouts get faint from 2 ~ 3 to 6 by a factor of l-7]!^03 t jBouwens et al.ll2006l ). On 
the other hand, UV-LFs of LAEs brighten by ~ 1 magnitude, if it comes from the pure 
luminosity evolution. This difference of LFs implies that Lya emitting population is more 
dominant at z = 6 than 2; = 3 at the given UV luminosity range. 



7.2. AGN fraction and Evolution of Faint AGN 

We identify 6 AGN-LAEs in total. In our photometric samples of LAEs (satisfying our 
color-selection criteria), we find that 3 out of 356 z = 3.1 LAEs and 1 out of 101 z = 3.7 
LAEs host AGNs (Section 14.61) . Thus, the fraction of AGNs is about 1% at z = 3.1 and 



3.7, which is consistent with the result of iGawiser et al.l (j2006[ ). Note that the detection 



limits of our LAEs reach (1 - 4) x lO''^ ergs s"^ (~ 0.2 - 0.4L^y^). This AGN fraction is the 
lower limit, since our multi- wavelength data are probably not deep enough to identify faint 
AGNs. We estimate the relation between Lya luminosity and the detection limits of our 
multi-wavelength data. By the comparison of our detection limits and the QSO templates in 
Figure UM (see Section 146) for the definition of QSO templates), our X-ray data likely reach 
the deepest detection limit for AGNs among our multi- wavelength data. We normalize the 
QSO templates to the detection limit of our X-ray data, and measure Lya luminosities of 
the QSO templates. We find Lya-luminosity limits of L[Lya) = 9.5 x 10^^ and 1.5 x 10^^ 
erg s~^ for 2; = 3.1 and 3.7 objects, respectively, assuming that the SED of AGN-LAEs is 
the same as the QSO templates. Down to these Lya-luminosity limits, we find 3 out of 25 
(1 out of 8) photometrically selected LAEs have an AGN at 2 = 3.1 (z = 3.7). Thus, the 
AGN fraction is estimated to be ~ 10% down to L(Lya) ~ 1 x 10^^ erg s~^ for z = 3.1 and 
3.7 LAEs. This Lya-luminosity limit corresponds to ^ 2L£y^. Note that the AGN fraction 
of LAEs is as small as 10% even with this bright detection limit of Lya luminosity. 

If we take a Lya limit for the very bright luminosity of logL(Lya) > 43.6 and > 43.4 
erg s~^ in our z = 3.1 and 3.7 samples, respectively, all of LAEs (2 out of 2 LAEs and 1 out 
of 1 LAE at 2; = 3.1 and 3.7) have AGN activities. These very bright limits correspond to 
the brightest bins of our LF in Figure [19] (see Section ISTT!) . There appear possible humps of 
LF at 2; = 3.1 and 3.7 which are produced by the AGN-LAEs. The number density in these 
brightest bins is about ~ 10~^ Mpc~^. Thus, a shallow wide-field LAE survey generally finds 
this AGN-LAE population. 

Steidel et al. (2002) have obtained the fraction of AGNs in their z = 3 LBG samples 
down to ^ M* + 1 with optical spectra. They have found the AGN fraction of 3% for their 
LBG sample with a relatively faint continuum limit of ^ M* + 1. Since the detection limits 
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of this study is different from ours, we cannot directly compare the AGN fraction of their 
LBGs and our LAEs. Since we find the AGN fraction (and its lower limit) of our LAEs 
is 10% (and > 1%) for our ~ 2L*y^ (and the entire [~ 0.2 - 0.4L*y^]) samples, the AGN 
fraction of LAEs and LBGs would be comparable. 

For 2; = 5.7 LAEs, we identify no AGN with our spectra or multi- wavelength data. In 
the same manner as the calculations above, we estimate a Lya detection limit of z = 5.7 
LAEs that corresponds to the X-ray luminosity limit with the QSO templates. We find 
that the Lya detection limit is L{Lya) ~ 4 x 10^^ erg s^^, but we have no LAEs with 
L{Lya) > 4 x 10^^ erg s~^. On the other hand, there are two AGN-LAEs that are brighter 
than L{Lya) ~ 4 x 10^^ erg s~^ (i.e. logL(LyQ;) > 43.6 erg s~^) at 2; = 3.1. Therefore, it 
may imply that LAEs with a bright AGN such found at z = 3 would disappear dX z = 5.7. In 
fact, no AGN signatures such as Nv emission have been confirmed in spectr a of the brightest 



2; = 5.7 LAEs that are found in other wide- field studies (jWestra et al.ll2006l ). No detection of 
AGN in 2 = 5.7 LAEs implies that the number density of LAEs with AGN activities would 
drop from 2; = 3 to 5.7. Recent QSO surveys have measured QSO luminosity functions 
accurately, and found that the space density of bright QSOs with Mttv < —26.7 decreases 



from 2: = 3 to 6 by more than an order of magnitude (e.g. ICroom et al.l |2004J : iFan et al. 



2OOII . I2OO4J ). This trend is the same as that of our finding in our LAE samples. Because 

our AGN-LAEs are as faint as Muv — —23 21, our AGN-LAEs probe for evolution of 

relatively faint AGN. Thus, evolution of faint AGNs appears to follow the evolution of QSOs. 
This faint-AGN evolution implies that not only bright QSOs but also faint AGNs can provide 
less significant contribution than star-forming galaxies for reionizing the Universe at z > 6. 
In other words, UV radiation from star formation likely becomes more dominant than that 
from faint AGNs at 2; ~ 6. 



7.3. Star-Formation Rates Indicated from Multi- Wavelength Data 

In Section HJ we have identified AGNs with our multi-wavelength data. The multi- 
wavelength data are too shallow to identify star-forming activities at high redshifts. The 
limiting SFRs of our multi-wavelength are as large as ~ lO^Mgyr"^ for our X-ray and 24/im 
data and ~ lOOOMoyr-^ for our radio data, respectively. However, we find that our radio 
data provide a meaningful limit on SFR of our LAEs by stacking analysis. We exclude LAEs 
with an individual radio det ection, and stack radio images of our photometrically selected 



LAEs in the same manner as lSimpson et al.l (l2006bl ). The stacked images show no detection 



in the radio band. The 3(T-upper limits of the stacked images are (/i.4ghz) = (1.8,7.0,2.7) 
/xJy for our z = (3.1, 3.7, 5.7) LAEs. 
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The relationship between SFR and radio luminosity (Li.4ghz) is given by 



SFR{MqYt- 



5.3 X 10^^{u/GHz) 



(11) 



for high- mass stars (M > 5Mq; ICondonlll992l : see also ISerjeant et al.ll2002l ). where u is the 
frequency and a is the non-thermal ratio spectral index {a ~ 0.8). Assuming this SFR v.s. 
radio-luminosity relation, we obtain 3a upper limits of (SFR) = (29, 162, 160) MqJt^^ for 
our z = (3.1,3.7,5.7) LAEs. Thus, the average SFRs of LAEs are smaller than ~ 30 and 
160 Mf^yr "-^ at z = 3.1 and 3. 7/5.7, respectively. The result of our z = 5.7 LAEs is similar 
to that of ICarilli et al.l (120061 ) who report the 2a upper limit of (/i.4ghz) = 2.5 fiJj (i.e. 
~ 100 M0yr~^) for their stacked radio image oi z = 5.7 LAEs in COSMOS field. The upper 
limits of radio SFR are consistent with the SFRs of typical (L*) LAEs, ^5 — 10 M0yr~^ 
(Section 15.31) . estimated from UV and Lya luminosities. This agreement between radio and 
optical SFRs implies that majority of LAEs are not associated with dusty starbursts with 
significant (xlO — 100) hidden star- formation. It is consistent with the fact that no LAEs 
are detected in submm bands (Section l4.4p . These results support the small-dust extinction 
of LAEs from the blue UV-continuum color (Section 17. 5p . On the other hand, no significant 
excess of radio SFR may reject the possibility of an extremely-top heavy initial mass function 
for star-formation of LAEs, since radio fluxes are originated from non-thermal radiation of 
supernovae at the final stage of massive stars. 

No radio detections in the stacked images also su ggest that typical L AEs do not harbor 



low-luminosity radio AGNs at z = 3.1 — 5.7 (see also ICarilli et al.l 120061 for z = 5.7 LAEs). 



The radio fluxes of our AGN-LAEs are /i.4ghz ^ 100 - 200/iJy (Table i]) at ;z = 3.1 and 3.7, 
while our 3a upper limits of stacked images are 1.8 — 7.0/iJy. Thus, typical LAEs do not 
have a low-luminosity radio AGN whose flux is brighter than 1/10 — 1/100 of typical radio 
AGN-LAEs at ^ = 3 - 4. 



7.4. Lya Emissivity 



Figure [25] shows that UV-bright (Muy < —21.5) galaxies do not have a large EWq 
{EWq > 80A) at 2; = 3 — 6. lAndo et al.l (120061 ) suggest that UV-bright galaxies are a dusty 



and chemically-evolved population, and that such UV-bright galaxies do not have a large EW. 
There are two other possibilities to explain this trend. First, the average stellar population 
of UV-bright galaxies would be older than that of UV-faint galaxies. Since Lya emission 
can be efficiently produced by hard-UV continuum from early-type stars, g alaxies can emi t 
a strong Lya emission of EWq ^ lOOA within ~ 10 Myr after a starburst ( Schaererl |2003| ) . 
Second, HI clouds of inter-stellar medium or inflows may be rich in UV bright population. 
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and absorb a Lya emission line. It is likely that the combination of age, neutral hydrogen, 
dust, and metallicity makes the deficit of large EW galaxies at bright UV magnitude. 



7.5. Implications from the Blue UV continuum 



We find that UV-continuum colors of LAEs are bluer than those of dropout galaxies in 
Section [631 This implies that LAEs are less-dusty and/or younger population than dropout 
galaxies. Although more detailed analysis with near-infrared data and population synthesis 
models will be presented in Ono e t al. (in prep a ration) , we compare the UV-continuum colors 



of LAEs with a simple model. lOuchi et al.l (l2004al ) have obtained the relation between 



i ' — z' color and dust ext inction, E{B — V\ with the stellar populatio n synthesis model 
( iBruzual fc Chariot I l2003l ) and the extinction law (ICalzetti et al.l 120001 ) which reproduce 
a typical UV-optical SED of dropout galaxies at z ~ 3. The relation of color and dust 
extinction is, 

E(5-V) ^0.0162 + L18(i'-^'), (12) 



Ouchi et al.l (j2004al ) also confirm that this relation agrees with the one of local starbursts. 
This relation assumes that the difference of UV color (or UV slope) is caused by dust ex- 
tinction, since the model fixes the age, metallicity, initial mass function, and star-formation 
history. We apply this simple model to our LAEs in Figure [221 Figure [211 implies that 
LAEs have an extinction of E{B — V)^ 0.05, while that dropout galaxies typically have 
E{B — V)^ 0.1 — 0.2. There are some very blue LAEs whose UV colors scatter around 
—0.2 — 0.0. The colors of these objects are still comparable to the color of E{B — V') = 0.0 
within a ^ 2cr level. However, if this very blue UV color is true, the colors of LAEs can be 
explained not only by e xtinction, but also by stellar population that is bluer than that of 
typical dropout galaxies. iGronwall et al.l (120071 ) estimate EiB — V^ of z = 3.1 LAEs with the 
independent approach, i.e. the comparison of UV and Lyct luminosity. They have obtained 
E{B — V) < 0.05, which is consistent with ours. 



8. Conclusions 

We have made a photometric sample of 858 LAE candidates and a spectroscopic sample 
of 84 confirmed LAEs at z = 3.1, 3.7, and 5.7 down to L(LyQ;) > 3 x 10"^^ erg s~^ in a 1 
deg^ sky of the Subaru XMM-Newton Deep Survey Field. Based on these samples, we have 
studied photometric and spectroscopic properties with the aid of multi- wavelength data from 
X-ray to radio. The major findings of our study are summarized below. 
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1. We derive the LF of Lya emission for each redshift and provide the best-fit Schechter 
parameters as summarized in Table 5, by carefully taking into account the statistical error 
and the field-to-field variation. We find that the apparent Lja LF shows no significant 
evolution between z = 3.1 and 5.7. Evolution, if any, is limited within factors of 1.8 and 
2.7 in L* and 0*, respectively. Although no evolution of the Lya LF over 3 < z < 5.7 has 
already been suggested by previous studies, we now have much stronger and more reliable 
constraints with smaller systematic and statistical errors. Since the average Lya opacity of 
the IGM is larger at earlier epochs, the absence of evolution in the apparent Lya LF implies 
an intrinsic brightening of Lya luminosity with increasing redshift being canceled out by 
increasing absorption of the IGM. We note, however, that this claim is true only if the two 
assumptions made in Section [5.1.21 are correct. 

2. The LF of UV continuum for LAEs shows an increase in number density and/or 
UV luminosity from 2; ~ 3 — 4 to 5.7. This evolution of the UV LF is consistent with 
the plausible evolution in the intrinsic Lya LF, since we find that the distribution of the 
rest-frame intrinsic EW of LAEs does not change or positively evolves with redshift. On 
the other hand, the UV LF of dropout galaxies is known to show an opposite evolutionary 
trend to that of LAEs, namely, a decrease in number density and/or UV luminosity with 
increasing redshift. This implies that the ratio in number density of LAEs to dropout galaxies 
increases from z ~ 3 — 4 to 5.7, and that galaxies with Lya emission are more common at 
earlier epochs. 

3. We identify six LAEs (in total) with AGN activities. In our photometric sample, 
three out of 356 z = 3.1 LAEs and one out of 101 z = 3.7 LAEs host AGNs. Thus, the 
lower limit of the AGN fraction down to (1 - 4) x 10"^^ g^g 8"^ (~ 0.2 - 0.4LJ^y^) is ~ 1% at 
z = 3.1 and 3.7. We use QSO templates given in the literature to find that ~ 10% of LAEs 
are AGNs above the Lya luminosity where our sample is complete, L(Lya) ~ 1 x 10^^ erg 
s-^ (~ 2LlyJ. It is notable that 100% of LAEs (2 out of 2 LAEs and 1 out of 1 LAE at 
z = 3.1 and 3.7) host AGNs at the bright end of the LF with logL(Lya) > 43.6 and > 43.4 
erg s~^ at z = 3.1 and 3.7, respectively. Shallow and wide-field narrow-band imaging would 
be an efficient way to search for LAEs hosting AGNs. 

4. We obtain the 3a upper limits of stacked radio images, (/i.4ghz) = (1-8, 7.0, 2.7) yuJy 
for our z = (3.1, 3.7, 5.7) LAEs. These upper limits indicate (SFR) = (29, 16 2 160) MqJt-^ 



at z = (3.1,3.7,5.7) based on the SFR-radio luminosity relation of ICondonI (119921 ). These 
upper limits are consistent with the SFRs estimated from Lya and UV luminosities. Thus, 
the majority of LAEs are not associated with dusty starbursts with significant (xlO — 100) 
hidden star-forming activities. No radio detections in the stacked images also indicate that 
typical LAEs do not harbor a low-luminosity radio AGN at z = 3.1—5.7. At z = 3—4, typical 
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LAEs do not have a low-luminosity radio AGN whose flux is brighter than 1/10 — 1/100 of 
typical radio AGN-LAEs of our SXDS field. 

5. The fraction of LAEs that have very large Lya EWs does not change significantly 
over 3.1 < 2 < 5.7; the fraction of objects with EW'^^ > 240A is 10-40 %. 

6. If we naively assume that all of the Lya or UV emission in our LAEs originate from 
star-forming activities, LAEs contribute about ~ 20 — 40% of the entire cosmic SFRDs at 
z = 3-6. 

7. We find a deficit of large-EW galaxies with bright UV cont inuum emission at a ll 
three redshi f ts. Th i s finding not only confir ms the previous results of IShapley et al.l (120031 ): 
Ando et al.l (120061 ): IShimasaku et al.l (120061 ) for z = 3 and 5 — 6, but also claims that this 
trend is common in the early universe. This deficit is probably caused by a combination of 
old ages and large amount of neutral hydrogen, dust, and metallicity in UV-bright galaxies. 

8. The UV-continuum colors of LAEs are generally bluer than those of dropout galaxies 
at the same UV luminosity. If we assume that LAEs have the same stellar population as 
dropout galaxies, then the UV spectral slope of LAEs indicates E{B — V) < 0.05. The 
very blue LAEs probably have extremely low dust extinction and/or very young stellar 
populations. 

9. We stack the spectra of the confirmed LAEs to estimate the 3a upper limit of 
fueu/ fhya to bc 0.02 for z = 3.1 and 0.06 for z = 3.7. These low fueii/ fhya values are 
consistent with no or little primordial population being included in our sample. 
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Table 1. Summary of Imaging Observations and Data 



Band 


Field Name(s) 


Exposure Time 


PSF 


size''' 


Area 




Date of Observations 






(sec) 


(arcsec) 


(arcmin^) 


(5(T AB mag) 






AfB503 


SXDS-C 




4200 


0.79 


(0.81) 


657 


25.3 


2003 


Oct 23 


NB503 


SXDS-N 




4200 


1.05 


(1.05) 


766 


25.1 


2003 


Oct 23 


NB503 


SXDS-S 




4200 


0.91 


(0.91) 


827 


25.3 


2003 


Oct 23 


NB503 


SXDS-E 




4200 


0.69 


(0.83) 


685 


25.4 


2003 


Oct 24 


NB503 


SXDS-W 




5400 


0.67 


(0.83) 


603 


25.5 


2003 


Oct 24 


NB570 


SXDS-C 




4200 


0.91 


(0.91) 


627 


24.6 


2003 


Oct 22,26 


NB570 


SXDS-N 




4200 


0.91 


(0.91) 


781 


24.9 


2003 


Oct 26 


NB570 


SXDS-S 




5400 


0.69 


(0.81) 


833 


25.2 


2003 


Oct 22 


NB570 


SXDS-E 




4200 


0.75 


(0.83) 


674 


24.9 


2003 


Oct 26 


NB570 


SXDS-W 




4200 


0.87 


(0.87) 


559 


24.8 


2003 


Oct 26 


NBSie 


SXDS-C 




17182 


0.65 


(0.81) 


676 


26.0 


2003 


Sep 28-30, Oct 22 


NB816 


SXDS-N 




14400 


0.73 


(0.85) 


810 


26.0 


2003 


Sep 29, Oct 22 


NBS16 


SXDS-S 




14400 


0.65 


(0.81) 


835 


26.1 


2003 


Sep 28,30 Oct 22 


NB816 


SXDS-E 




14255 


0.67 


(0.83) 


722 


26.0 


2003 


Sep 29-30, Oct 22 


NB816 


SXDS-W 




20400 


0.69 


(0.83) 


678 


25.9 


2003 


Sep 30, Oct 22 


Archival broad-band data^t. 


B 


SXDS-C,N,S,E,W 


19800 


- 20700 


0.78- 


-0.84 


915 - 979 


27.5 - 27.8 






V 


SXDS-C,N,S,E,W 


17460 


- 19260 


0.72 - 


- 0.82 


915 - 979 


27.1 - 27.2 






R 


SXDS-C,N,S,E,W 


13920 


- 14880 


0.74- 


- 0.82 


915 - 979 


27.0 - 27.2 






i' 


SXDS-C,N,S,E,W 


18540 


- 38820 


0.68 - 


- 0.82 


915 - 979 


26.9 - 27.1 






z' 


SXDS-C,N,S,E,W 


11040 


- 18660 


0.70 - 


- 0.76 


915 - 979 


25.8 - 26.1 







-t-Thc FWHM of PSFs in the reduced image. The value in parenthesis indicates the FWHM of PSF that is matched 
with broad-band images in each field. 

tThe limiting magnitude defined by a 5<t sky noise in a 2"-diameter circular aperture. 

ttThe archival broad-band data of SXDS presented in Furusawa et al. (in preparation). We summarize the properties 
of the 5-field images on a single line. Note that the exposure time is not a total of the 5 fields, but 1 field. More details 
are presented in Table 2 of Furusawa et al. in preparation. 



Table 2. Lya Emitters with spectroscopic redshifts 



Object Name a(J2000)" <5(J2000)" tonb msB Snb^ Sbb^ L(Lya) EW^'^'^ Muv SFR Note 

(mag) (mag) (kpc) (kpc) (lO^^erg s"!) (A) (mag) (Moyr"!) 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 



: 3.1 LAEs 



NB503-C-112924 


02 


18 


50 


-04 


50 


54 


3 


134 


24 


9 


27 


8 


4.1 




3.0 ±0.5 


264.6lg4^3^ 


-18 





0.9 


NB503-C-114646 


02 


17 


44 


-04 


50 


32 


3 


135 


23 


9 


25 


4 


4.1 


<4 


6.6 ±0.5 


55.9^g 


-20 


2 


6.5 


NB503-C-15140 


02 


18 


58 


-05 


12 


26 


3 


147 


23 


6 


25 


6 


4.4 


<4 


9.0 ±0.8 


62.9+11 


-20 


2 


6.5 


NB503-C-29891 


02 


17 


21 


-05 


09 


12 


3 


105 


23 


7 


25 


9 


7.1 


<4 


14.8 ±2.1 


124.8ii?-^ 


-19 


7 


4.2 


NB503-C-33149 


02 


18 


39 


-05 


08 


27 


3 


141 


24 


6 


26 


6 


6.9 


4.5 


3.5 ±0.6 




-19 





2.1 


NB503-C-46996 


02 


18 


33 


-05 


04 


57 


3 


129 


23 


9 


25 


6 


4.3 


4.3 


5.9 ±0.7 




-20 





5.4 


NB503-C-49497 


02 


18 


06 


-05 


04 


08 


3 


129 


21 


8 


24 


4 


8.3 


6.2 


44.4 ± 1.3 


62.6+ig 


-21 


7 


(25.3) 


NB503-C-51702 


02 


18 


07 


-05 


03 


49 


3 


125 


24 


4 


25 


8 


5.4 


<4 


3.5 ±0.7 




-20 


1 


6.2 


NB503-C-52511 


02 


17 


12 


-05 


03 


37 


3 


125 


24 


7 


26 


8 


5.6 


<4 


3.2 ±0.5 


104.3t^j5;| 


-18 


8 


1.7 


NB503-C-53382 


02 


17 


05 


-05 


03 


25 


3. 


152: 


24 


5 


25 


5 


4.7 


4.6 


4.4 ±0.5 


4i.8ig:5 


-20 


1 


6.2 


NB503-C-57309 


02 


17 


08 


-05 


02 


26 


3 


136 


24 


6 


25 


8 


7.3 


4.3 


3.3 ±0.5 




-20 


1 


6.0 


NB503-C-66311 


02 


17 


27 


-05 


00 


21 


3 


144 


24 


6 


26 


4 


5.1 


<4 


3.5 ±0.5 


69.1+1^-° 


-19 


2 


2.7 


NB503-C-88741 


02 


17 


00 


-04 


55 


40 


3 


113 


23 


6 


25 


6 


<4 


<4 


11.7± 0.7 


128.2tlJ-* 


-20 


1 


6.2 


NB503-C-90464 


02 


17 


43 


-04 


55 


21 


3 


111 


24 


8 


27 


1 


4.1 




4.2 ±0.8 


204.2^53 g 


-18 


6 


1.6 


NB503-C-99469 


02 


18 


03 


-04 


53 


31 


3 


133 


24 


9 


27 


2 


4.5 




2.7 ±0.5 


137.2+^:^"^ 


-18 


4 


1.2 


NB503-E-161419 


02 


20 


24 


-04 


52 


05 


3 


105 


24 


6 


26 


3 


4.6 


4.2 


6.1 ± 1.8 


86.2+2|0 


-19 


4 


3.2 


NB503-E-49946 


02 


19 


09 


-04 


47 


22 


3 


152 


24 


3 


26 





4.7 


4.6 


5.2 ±0.7 


65.2+11 


-19 


5 


3.5 


NB503-E-60750 


02 


19 


16 


-04 


58 


38 


3 


129 


23 


5 


25 


1 


<4 


<4 


9.2 ±0.6 


58-4+39 


-20 


5 


8.4 


NB503-N-32683 


02 


17 


41 


-04 


43 


47 


3 


106 


24 


4 


25 


9 


5.0 


7.4 


7.3 ± 1.0 


114.1+in 


-19 


8 


4.6 


NB503-N-33460 


02 


17 


27 


-04 


43 


36 


3 


119 


23 


3 


25 


6 


5.3 


4.1 


12.9 ±0.6 


121.3l«:^ 


-20 


1 


6.0 


NB503-N-34067 


02 


17 


44 


-04 


43 


29 


3 


152 


24 


6 


> 28.0 


<4 




4.4 ±0.7 


'iA'i r>-t-292.7 


-17 


6 


0.6 


NB503-N-35820 


02 


17 


35 


-04 


42 


59 


3 


102 


23 


6 


25 


3 


6.6 


6.1 


21.1 ± 1.9 


1 c;+12.8 
loo.o_j2.8 


-20 


8 


(11.0) 


NB503-N-42377 


02 


17 


42 


-04 


41 


08 


3 


154 


22 


2 


23 


9 


4.2 


4.4 


38.8 ±0.8 


73.6+11 


-21 


8 


29.0 


NB503-N-72492 


02 


17 


42 


-04 


33 


20 


3 


108 


24 


3 


26 


2 


7.2 


6.3 


8.0 ± 1.0 


146.8ii:^ 


-19 


4 


3.2 


NB503-N-80475 


02 


17 


41 


-04 


31 


30 


3 


123 


21 


8 


22 


6 


<4 


4.3 


31.0 ±0.6 


lS.6+J,i 


-23 


2 


(100.7) 


NB503-N-87126 


02 


17 


54 


-04 


30 


06 


3 


112 


23 


5 


25 


4 


6.1 


4.6 


12.4 ±1.1 


83.4+^1 


-20 


3 


7.0 


NB503-N-92276 


02 


17 


42 


-04 


28 


58 


3 


133 


24 


3 


26 


4 


4.2 


5.6 


4.5 ±0.7 


7r q+14.7 
'^•^-11.9 


-19 


7 


4.0 


NB503-S-105705 


02 


17 


06 


-05 


16 


34 


3 


106 


24 


6 


26 


5 


7.1 


5.3 


6.6 ± 1.3 




-19 


3 


2.8 


NB503-S-35558 


02 


18 


18 


-05 


31 


57 


3 


132 


24 


7 


26 


6 


6.6 


5.4 


3.4 ±0.5 




-19 


6 


3.7 


NB503-S-4r3244 


02 


18 


26 


-05 


29 


45 


3 


156 


23 


1 


24 


8 


<4 


<4 


18.5 ± 0.6 


82.3+^-^ 


-21 


1 


15.7 


NB503-S-54416 


02 


17 


05 


-05 


27 


35 


3 


133 


23 


9 


26 





5.9 


9.3 


6.7 ±0.6 




-19 


6 


3.7 


NB503-S-56809 


02 


17 


57 


-05 


27 


00 


3 


123 


24 


2 


27 


3 


5.9 




5.9 ±0.8 


233.7+§^-? 


-18 


3 


1.2 


NB503-S-66012 


02 


17 


23 


-05 


24 


49 


3 


167 


24 


6 


25 


6 


4.3 


6.4 


7.8 ± 1.0 


73 4+10.5 


-20 





5.6 


NB503-S-89393 


02 


19 


04 


-05 


19 


58 


3 


125 


24 


4 


26 





5.9 


5.9 


3.8 ±0.8 


42.6tr2 


-19 


6 


3.8 


NB503-S-94275 


02 


18 


41 


-05 


18 


49 


3 


102 


22 


8 


24 


5 


4.1 


<4 


46.5 ± 1.9 


164.618:3 


-21 


2 


16.1 


NB503-W-100818 


02 


16 


04 


-04 


49 


43 


3 


125 


24 





26 





4.1 


<4 


6.3 ±0.6 


85.0llJi 


-19 


6 


3.8 


NB503-W-31855 


02 


16 


52 


-04 


56 


43 


3 


130 


24 


4 


26 


9 


5.1 


9.1 


4.6 ±0.6 


1 Q-t-39.2 
lOI .o_27^g 


-18 


6 


1.5 


NB503-W-67066 


02 


16 


27 


-05 


02 


43 


3 


139 


23 


9 


25 


8 


7.3 


11.9 


6.3 ±0.8 




-19 


6 


3.7 


NB503-W-85295 


02 


16 


14 


-05 


02 


15 


3 


141 


24 


5 


26 





<4 


4.3 


3.8 ±0.6 


51-2+8.2 


-19 


6 


3.7 


NB503-W-97848 


02 


16 


06 


-04 


48 


12 


3 


121 


24 


6 


26 


6 


4.1 


<4 


3.7 ±0.6 


ini 1+21-5 


-19 


1 


2.3 


NB503-W-98167 


02 


16 


06 


-04 


50 


57 


3 


124 


23 


8 


25 


6 


4.1 


<4 


6.9 ±0.7 


63.5l^J 


-20 





5.4 



AGN 



AGN 
Civ?!' 

AGN 
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Table 2 — Continued 



Object Name 
(1) 



a(J2000)" (5(J2000)" 2* 



(2) 



(3) 



(4) 



m-NB 
(mag) 
(5) 



niBB 
(mag) 
(6) 



Snb'' Sbb'' L(Lya) 
(kpc) (kpc) (lO^'^erg s-i) 
(7) (8) (9) 



(A) 
(10) 



Mvv SFR 
(mag) (Moyr-l) 
(11) (12) 



Note 



(13) 



z = 3.7 LAEs 



NB570-C-in7443 


02 


17 


47 


-04 


52 


30 


3 


648 


23 


4 


24 


3 


4.7 


< 4 


33.6 ± 4.7 


75 2+i°o* 


-21 


7 


26.7 


NB570-C-1 16152 


02 


18 


56 


-04 


50 


59 


3 


683 


24 


5 


26 


4 


4.1 


< 4 


4.7 ± 0.8 


87 3+JI'2 
°' -^-16.5 


-19 


g 


3.7 


NB57n-C-l 3685 


02 


17 


09 


-05 


12 


29 


3 


691 


23 


9 


27 


Q 


5.7 




8.3 ±1.1 


335 Q+i^?;? 


-18 


3 


1.1 


NB570-C-30225 


02 


18 


26 


-05 


10 


03 


3 


699 


23 





24 


8 


5.7 


4.7 


19.0 ± 1.1 


57 7+o'S 


-21 


2 


16.1 


NB570-C-40856 


02 


17 


01 


-05 


07 


28 


3 


686 


24 





26 


8 


4.8 


4.5 


7.2 ± 0.9 


161.8+?Si 

— o\j.Z 


-19 


2 


2.6 


NB570-C- 78067 


02 


17 


13 


-04 


58 


53 


3 


684 


23 


5 


24 


6 


5.7 


4.5 


7.9 ± 1.6 




-21 


1 


15.5 


NB570-E-103151 


02 


19 


44 


-04 


50 


55 


3 


669 


23 


2 


25 


8 


4.6 


<4 


21.9 ±2.0 


150.0liJ:J 


-20 


1 


6.0 


NB570-E-142687 


02 


20 


11 


-04 


50 


08 


3 


644 


24 


4 


25 


4 


5.6 


4.2 


16.6 ± 7.0 


95.4if,:? 


-20 


5 


8.8 


NB570-E-169608 


02 


20 


30 


-04 


48 


15 


3 


702 


24 


6 


26 


8 


5.3 


5.0 


4.7 ±0.9 


110.412^:6 


-19 


2 


2.7 


NB570-E-58808 


02 


19 


13 


-04 


51 


59 


3 


692 


23 


3 


25 


7 


<4 


<4 


14.2 ±0.9 


1 26 4+111 


-20 


4 


7.8 


NB570-E-62593 


02 


19 


15 


-04 


55 


11 


3 


672 


23 


6 


25 


5 


6.0 


<4 


12.8 ± 1.7 


67.41^:2 


-20 


7 


10.6 


NB570-E-65332 


02 


19 


17 


-05 


07 


39 


3 


693 


23 


8 


25 


9 


<4 


<4 


9.1 ±0.8 


92.8li°5^ 


-20 


1 


6.1 


NB570-N-33608 


02 


17 


27 


-04 


44 


14 


3 


683 


24 





25 


9 


5.7 


<4 


7.1 ± 1.1 


63.1+iV 


-20 





5.5 


NB570-N-34841 


02 


18 


31 


-04 


43 


54 


3 


724 


21 


8 


22 


2 


<4 


<4 


145.0 ± 1.5 


48.9+«'? 


-23 


8 


(178.2) 


NB570-N-42463 


02 


18 


40 


-04 


42 


11 


3 


692 


24 


4 


26 


4 


4.4 


<4 


4.7 ±0.8 


82.7111-1 


-19 


6 


3.9 


NB570-N-86993 


02 


18 


13 


-04 


30 


57 


3 


668 


23 


8 


25 


7 


5.2 


<4 


11.5 ±1.5 


90.7112.8 


-20 


4 


7.6 


NB570-S-125887 


02 


18 


07 


-05 


20 


48 


3 


638 


24 


3 


25 


2 


5.9 


4.4 


28.4 ± 11.4 




-20 


8 


11.0 


NB570-S-132143 


02 


18 


27 


-05 


19 


47 


3 


677 


23 


4 


25 


9 


4.1 


<4 


1-1.2 ± 1.3 


ii8.9ll*:i 


-20 





5.5 


NB570-S-141871 


02 


18 


10 


-05 


18 


13 


3 


684 


24 


5 


26 


9 


<4 


<4 


4.6 ±0.8 


134.4lg:^ 


-18 


9 


2.1 


NB570-S- 154087 


02 


17 


43 


-05 


16 


11 


3 


682 


24 


9 


27 


3 






3.3 ±0.9 


141 Q+72.6 


-18 


4 


1.2 


NB570-S-175068 


02 


18 


28 


-05 


13 


06 


3 


701 


24 


8 


26 


8 


5.7 


<4 


3.8 ± 1.0 


84.1126-8 


-19 


1 


2.3 


NB570-S-32766 


02 


18 


35 


-05 


35 


50 


3 


671 


23 





25 


4 


<4 


<4 


24.8 ±1.5 


129.711:1 


-20 


7 


10.6 


NB570-S-84321 


02 


17 


45 


-05 


27 


35 


3 


648 


24 


4 


25 


1 


7.2 


<4 


12.4 ± 4.6 


01- ' -22.0 


-20 


9 


11.9 


NB570-S-99194 


02 


17 


14 


-05 


25 


16 


3 


729 


24 


9 


25 


6 




<4 


11.6 ±2.7 


94.5123:1 


-20 


4 


7.9 


NB570-W-53415 


02 


16 


42 


-04 


58 


55 


3 


667 


24 


2 


25 


7 


<4 


<4 


8.2 ± 1.5 


66.3lli:« 


-20 


5 


8.4 


NB570-W-55371 


02 


16 


40 


-05 


01 


29 


3 


699 


24 


2 


24 


7 


4.7 


<4 


5.2 ±0.8 


i9.9li:°3 


-21 


2 


17.0 



AGN 



: 5.7 LAEs 



NB816-E-127266 


02 


20 


12 


-04 


49 


50 


5.681 


24 


9 


> 27.0 


<3 




8.0 ±1.3 


159.6+222^1 


-19 


2 


2.6 


NB816-E-129103 


02 


20 


13 


-04 


51 


09 


5.744 


24 


8 


26.1 


<3 




11.5 ±1.4 


64.3l«;i 


-20 


5 


8.8 


NB816-E-141288 


02 


20 


21 


-04 


53 


14 


5.671 


24 


8 


> 27.0 


<3 




9.6 ± 1.8 


120.6+42'ig3 


-19 


6 


3.9 


NB816-E-147538 


02 


20 


26 


-04 


52 


34 


5.718 


23 


7 


25.1 


conf. 


conf. 


18.9 ± 1.7 


29.2112 


-21 


4 


(20.1) 


NB816-S-36496 


02 


18 


22 


-05 


33 


37 


5.650 


25 


1 


26.5 


<3 




13.2 ±3.4 


133.2+^2.^ 


-20 


1 


6.1 


NB816-S-39206 


02 


18 


19 


-05 


33 


11 


5.676 


25 


3 


> 27.0 


4.5 




5.7 ± 1.5 


lot: 1+300.8 

10d.l_gg g 


-18 


9 


2.0 


NB816-S-41408 


02 


18 


14 


-05 


32 


49 


5.673 


23 


9 


25.7 


conf. 


conf. 


17.5 ±3.0 


42.71^2^8 


-21 


3 


(18.5) 


NB816-S-44568 


02 


18 


17 


-05 


32 


22 


5.644 


24 


6 


> 27.0 


<3 




36.0 ±4.7 


ocrc: f;+1327.0 
odO.O_4Q2 Q 


-19 





2.2 


NB816-S-46148 


02 


18 


23 


-05 


32 


05 


5.680 


24 


9 


> 27.0 


3.6 




8.2 ± 1.5 


1 74 0+311.2 


-18 


9 


2.0 


NB816-S-49611 


02 


17 


43 


-05 


31 


35 


5.629 


25 





25.7 


<3 


<3 


28.9 ± 14.2 


i26.il«^:! 


-21 


1 


15.7 


NB816-S-50308 


02 


17 


48 


-05 


31 


27 


5.690 


24 





26.0 


3.2 




15.0 ± 1.2 


81.41?^:^ 

91 0+9.5 
^^•■^-6.2 


-20 


7 


10.2 


NB816-S-59282 


02 


17 


51 


-05 


30 


03 


5.712: 


25 


1 


26.0 


3.8 




4.8 ± 1.3 


-20 


8 


10.8 
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Table 2 — Continued 



Object Name 
(1) 


a(J2000)'' 
(2) 


(5(J2000)'' 
(3) 


(4) 


(mag) 
(5) 


JTtBB 
(mag) 
(6) 


•S'nb^ 
(kpc) 
(7) 


(kpc) 
(8) 


L(LyQ) 
(lO^^erg s-i) 
(9) 


(A) 
(10) 


Muv 
(mag) 

(11) 


SFR Note 
(A%yr-i) 

(12) (13) 


NB816-S-61269 


02:17:45 


-05:29:36 


5.688 


23.7 


26.1 


3.5 




21.3 ± 1.7 


93.3l?^j 


-21.4 


19.2 


NB816-S-66352 


02:17:49 


-05:28:54 


5.696 


24.0 


25.9 


<3 


conf. 


16.2 ± 1.1 




-21.0 


(13.4) 


NB816-S-67673 


02:17:45 


-05:28:42 


5.751 


25.0 


26.1 


3.9 




10.9 ± 1.4 


61.6121? 


-20.4 


7.6 


NB816-S-70769 


02:17:43 


-05:28:07 


5.685: 


23.9 


26.1 


<3 




17.8 ± 1.2 


122.9l«;« 


-20.5 


8.5 


NB816-S-77389 


02:17:50 


-05:27:08 


5.693 


24.3 


> 27.0 


6.1 




11.9 ± 1.9 


106 4+107.0 


-19.6 


3.6 



Note. — (1): Object name. (2)-(3): RA and Dec. (4): Redshift. (5)-(6): Magnitudes in narrow and broad bands, 
(ArB503,i?), (NB570,R), and (AfB816,^'). for ^ = 3.1, 3.7, and 5.7 LAEs, respectively. (7)-(8) Half of FWHMs in physical kpc 
in the narrow and broad bands. We only show reliable FWHMs for sources with a > 5o" significance. (9): Lya luminosity in 
10^2 erg s~i. (10): Rest-frame apparent equivalent width of Lya emission line. (11): UV total magnitude defined in Section 
15.21 (12): Star-formation rate estimated from the UV total magnitude. The values in parenthesis present formal SFRs which 
are estimated for AGNs or confused objects. (13): AGN classification and notes. 

^See the published version for the exact coordinates with decimals. 

tin column (4), ":" marks objects with an uncertain spectroscopic identification. 

t "conf" in the columns (7) and (8) indicates the sources confused by (a) neighboring objoct(s), which do not give an accurate 
FWHM. 

^The object with a marginal Civ detection, NB503-N-42377. See Section liTTI 
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Table 3. Samples of Lya Emitters 



Redshift 


Survey Area 
(arcmin^) 


Magnitude Range! 
(AB mag) 


L, t 
(erg s-l) 


(A) 




-^spcc 


Selection Criteria 


2 = 3.1 ± 0.03 


3538 


Ar_B503 = 22.4 - 25.3 


> 1 X 10*2 


> 64 


356 


41 


eq.Q 


z = 3.7 ±0.03 


3474 


Ar_B570 = 22.7 - 24.7 


> 4 X 10''^ 


> 44 


101 


26 


eq.© 


z = 5.7 ±0.05 


3722 


Ar_B816 = 23.6 - 26.0 


> 3 X 10*2 


> 27 


401 


17 


eq.® 



■'■2"-diameter aperture magnitudes. 

tThe approximate limit of Lya luminosity. 

ttThe approximate limit of rest-frame apparent equivalent width of Lyo emission. The corresponding intrinsic equivalent 
widths corrected for the IGM absorption are estimated to be EW^"^^ = 80, 60, and 50A, for z = 3.1, 3.7, and 5.7 samples, 
respectively. 



Table 4. Lya Emitters with AGN activities 



ID(opt) RA Dec z Dot ID(X) sGp(X) 


/2-iOkoV 


sep(IR) 






ID(R) 


sep(R) /i.4GHz 




Muv 


FWHM(NB) Type 


{J2000) {J2000) 


{") 


(10~i5(,gs) 


{") 


(/^Jy) 


(mjy) 




(") 




(10*3 erg s- 


1) (mag) 


(kpc) 




(1) (2) (3) (4) (5) (6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


(18) 


NB503-C-4949702:18;06.826-05;04:08.803.129SPXooo X0700 


2.6 


7.1 ± 0.3 




< 240 


<4 






< 60 


4.4 


-21.67 


16.5 


BL 


NB503-N-3582002:17:35. 593-04:42:59.833. 102 SPoIoo ■ ■ ■ 




< 0.7 


0.9 


463 ± 48 








< 60 


2.1 


-20.77 


13.2 


BL? 


NB503-N-5538002:17:01. 676-04:37:20.24 (3.1) oPoIoR ■ ■ ■ 




< 0.6 


1.6 


378 ± 48 




VLA0267 


1.8 


190 ± 20 


4.2 


-20.64 


21.9 




NB503-N-8047502:17:41. 987-04:31:30.553. 123 SooIoR • • ■ 




< 0.9 


1.2 


431 ± 48 




VLA0351 


0.7 


147 ± 20 


3.1 


-23.17 


7.8 


BL 


NB570-N-7184202:18:57.776-04:34:31.53 (3.7) oPXooo X0950 


1.7 


4.4 ± 1.4 




< 240 








< 60 


3.0 


-21.65 


8.5 




NB570-N-3484102:18:31. 379-04:43:54.773. 724 Sooooo ■ ■ ■ 




< 1.3 




< 240 








< 60 


14.5 


-23.79 


< 6.1 


BL 



Note. — (1): Object name from the optical data. (2)-(3): RA and Dec in the optical data. (4): Redshift. The values in parenthesis are photometric redshifts. (5) The 
identification of AGN-LAEs in our spectroscopy and multi- wavelength data. "X", "I", "M", or "R" means that the LAE has either an X-ray, infrared, submm, or radio 
detection. The flags of "S" and "P" indicate a spectroscopically identified AGN and an object satisfying our color-selection criteria of LAEs, respectively, ("o" means no 
detection or identification in the corresponding condition.) (6) and (12): ID names in the X-ray and radio data presented in Ueda et al. (in preparation) and lSimpson et al.l 
l l2006al ). respectively. (7), (9), and (13): Separation between positions of the optical center and a counterpart of X-ray, infrared, and radio. (8), (10), (11), and (14): Fluxes of 
X-ray 2-10 keV (in 10~^^ erg s"'^ cm"'^), infrared 24/.tm (in ^tJy), submm 850/im (in mJy), and radio 1.4GHz (in /iJy), respectively. For the upper limits of X-ray fiuxes, we 
assume an absorption of neutral hydrogen with a column density of lO^'^ cm"'^ and a spectral power of 7 = 1.8. We omit the effects of X-ray reflection, which provides the 
most conservative upper limit in column (8). (15) Lya luminosity in 10*^ erg s~^ derived from our photometric data. (16) Total UV magnitude in AB mag. (17) FWHM in 
the narrow band in physical kpc. (18) Type of object. "BL" indicates an AGN with a broad-line emission (fpwHM ^ 1000 km s"'^) in high ionized lines (i.e. Nv, Siiv, Civ, 
Hell, and/or Cm]). NB503-N-35820 shows a possible broad line with 1321 ± 362 km s~^ in Civ, although the signal-to-noisc ratio is as small as 4. 
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Table 5. Summary of the Lya luminosity functions. 



(10-*Mpc- 


3)(10*2crg s- 


a o"EW Xr 
1) (lO^A) 1 


^obs 

;io-''Mpc- 


obs 

3)(1039erg s-^MpC^ 


tot 
f^Lya 

-3)(1039erg s-lJvIpc- 


</'o 

-3)(10-4Mpc- 


/'OLya 
=*)(1039erg s-lMpc-3) 


(1) 


(2) 


(3) 


(4) (5) (6) 


(7) 


(8) 


(9) 


(10) 


(11) 




a = -1.5 (fix) 


3.1 
3.7 
5.7 


9.2^2-5 
o 4+1.0 

7 7+'^-* 
'•'-3.9 


c- O+0.9 
0-°-0.7 

10.2t^« 


-1.5 1.3to ;Jl.57 
-1.5 1.5to40.70 
-1.5 2.7(fix)1.04 


15.0t^;2 
2 q+0.8 




9.4+2-0 

9.2tl:? 


io.8«:« 
3.7ll:S 

7 7+7.5 
'■ '-3.9 


ll-'J-2.0 
"•'-1.4 
^•^-3.7 


a = -1.0 (fix) 


3.1 
3.7 
5.7 


14.9tl^ 
c 7+2.6 


7 9+1.3 
'■^-1.6 

5.4irr 


-1.0 1.3toil.72 
-1.0 

-1.0 2.7(fix)1.10 




4 5+1-^ 
^■0-0. 8 
2 4+0.9 

q c:+2.3 


''•1-1.0 

4 1+1.3 

5 fi+3-6 
3-"-1.6 


17.5+I-? 


7.llli 

c 7+3.6 
^•'-1.6 


a = -2.0 (fix) 


3.1 
3.7 
5.7 


1 n+0.5 

■^•"-2.5 


9 1+2.6 
16.2t|i 


-2.0 1.3^0 il. 51 
-2.0 1.5^030.71 
-2.0 2.7(fix)0.97 


18.5«;? 
2.9«J 
o.a_4.0 


5 4+2-1 

r, O + 0.8 

^■■^-0.5 
o 4+4.8 




4.511:? 

-1 c+0.6 
1-3-0.3 
q c+3.5 
''-0-2.5 





Note. — (1): Redshift. (2)-(4): Best-fit Schechter parameters in units of 10 '^Mpc ^ and lO'^^erg s i for (p* and ^Lya' respectively. 
a is fixed to —1.5, —1.0, and —2.0 in 1-3, 4-6, 7-9 lines, respectively. (5): Best-fit iTew in units of 10^ A. For z = 5.7 LAEs, we fix 
"■ew = 270A., (6): Reduced of the fitting. (7)-(8): Number densities (in 10~^Mpc~3) and Lya luminosity densities (in lO'^^erg s~i 
Mpc"-^) calculated with the best-fit Schechter parameters down to the observed limit of Lya luminosity, i.e. logLxjya = 42.1, 42.6, and 
42.4 (erg s~i), for z = 3.1, 3.7, 5,7, respectively. (9): Inferred total Lya luminosity densities integrated down to LLya = with the 
best-fit Schechter parameters. (10): Same as (2) but for all LAEs with EW > OA estimated from the simulations. (11): Same as (9) but 
for all LAEs with EW > OA estimated from the simulations. 
We refer to the results of a = —1.5 as the best estimates of LPs. 

^ ^Lya apparent value, i.e. observed Lya luminosity with no correction for IGM absorption. If wc assume the simple IGM 

absorption model in Section 15.1.21 the apparent L^^^ is divided by a factor of (0.81, 0.73, 0.54) to obtain the intrinsic ^Lya at 2; = 
(3.1, 3.7, 5.7). For examples, in the case of a = —1.5, the intrinsic L^y^ is estimated to be (7.1+o'g, 14.1+2'ji 12.6+| g) X lO'^^erg s~i at 
z = (3.1,3.7,5.7). 
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Table 6. Summary of the UV luminosity functions. 



z 

(1) 


r 

(10-*Mpc-3) 
(2) 


M^500 
(mag) 

(3) 


(4) (5) 


Mag. Range 
(mag) 
(6) 


„obs 

(10-4Mpc-3) (lO^Serg 
(7) 


s-i Hz-iMpc-3) (lO^^erg 
(8) 


s-i Hz-iMpc-3) 
(9) 


3.1 
3.7 
5.7 


K 9+0.8 


-19.8 ±0.4 
-19.8(fix)t 
-20.6 ±0.6 


-1.6 0.43 
-1.6 2.74 
-1.6 0.29 


-21.9 < M < -18.9 
-21.7 < M < -19.7 
-22.1 < M < -20.6 


4.0t|I 


1 2+1-'' 

"i+o -i 
•J-^-o.i 

0.9tg;f 


4.4l^i 

4.it°o:i 



Note. — (1): Rcdshift. (2)-(4): Best-fit Schcchtcr parameters in units of lO^^Mpc^'' and AB magnitude for 4>* a-nd -'^^isoo- 
fixed to —1.6. (5): Reduced of the fitting. (6): Magnitude range of UV LFs that are used for the fitting. (7)- (8): Number densities 
(in 10~*Mpc~^) and UV luminosity densities (in lO^^erg Hz~^ Mpc~^) calculated with the best-fit Schechter parameters down 
to the observed limit of UV luminosity, i.e. M1500 = —18.9, —19.7, and —20.6, for z = 3.1, 3.7, 5.7, respectively. (9); Inferred total 
UV luminosity densities integrated down to M1500 = 00 with the best-fit Schechter parameters. 

tPor the z = 3.7 LF, we fix A^*5Qo to -19.8 that is the best estimate of our z = 3.1 LF. 
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Fig. 1. — Top panel : Response of the narrow-band filters, NB503, NB570, and NB816 (black 
solid lines), of this work, together with that of bro ad-band filters of B,VM4\ and z' (dotted 
lines) and two other narrow-band filters of NB711 (jOuchi et al.ll2003l : IShimasaku et al.ll2003l ) 
and NB921 (gray solid lines). A^i?921 filter is used for our on-going observations in the 
SXDS field. These response curves include atmospheric absorption, quantum efficiency, 
and transmittance of optical elements of the telescope and instrument. Bottom panel : The 
spectrum of night-sky emission. Note that the passbands of the narrow-bands do not include 
a strong sky emission. 
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Fig. 2. — Redshift distribution of our objects with a spectrum. Thin sohd hues in all three 
panels show the distribution of our LAEs and objects of the SXDS spectroscopic catalog. 
The shaded histograms indicate the sources detected at a 5a level with A^i?503 (top), NB570 
(middle), and NB816 (bottom). The thick solid lines in top, middle, and bottom panels 
represent our spectroscopic LAEs that meet the color criteria of equations ([I]), ([2]), and ([3]), 
respectively. 
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Fig. 3. — Spectra and snapshots of our z = 3.1 LAEs. Each object has a spectrum in right 
and snapshots of B, V, R, i', z', and A^i?503 images in left. Each snapshot is presented 
in a 6" x 6" box. In the panel of spectrum, the tick of y axis is marked in 2.5 x 10~^^ erg 
s~^ cm~^ for panels in which a factor is shown in the upper left corner, multiply the 

scale by this factor to obtain a correct scale. The wavelength (in x axis) is in unit of A. The 
object name and redshift (+AGN classification, if any) are presented in the right corner of 
each spectrum panel. 'P' ('N') in parentheses indicates that the object is (or is not) selected 
with the color criteria and included in our photometric sample (Section 13. ip . The vertical 
dotted lines mean the center of the emission line. The right bottom panel shows a typical 
spectrum of the sky background with an arbitrary normalization. 
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Fig. 4. — Same as Figure [HI but for our z = 3.7 LAEs. 
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Fig. 5. — Same as Figure [3], but for our z = 5.7 LAEs. The ticks of y axis in a spectrum 
panel are marked in 0.5 x 10~^^ erg s~^ cm~^ A"^. 
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Fig. 6. — The redshift distribution of spectroscopically identified LAEs at z = 3.1 (left), 
3.7(middle), and 5.7 (right). Thin sohd fines indicate aU tfie spectroscopicaUy identified 
objects. The shaded histograms present the sources at the 5a level in NB503 (left), NB570 
(middle), and NB81Q (right) bands. The thick solid lines are histograms for LAEs selected 
by the color criteria of equations ([1]) (left), ([2]) (middle), and ([3]) (right). The dashed lines 
represent selection functions of LAEs that are simply calculated from the response curves 
of narrow-band filters. The dotted lines are selection functions from the results of our 
simulations in Section 15. 1[ Both of the selection functions are normalized by the numbers of 
identified LAEs. 
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Fig. 7. — Color magnitude diagrams of A^i?503 (left), NB570 (middle), and NB816 (right). 
The black dots plot colors of all the detected objects. The photometrically-selected LAEs 
are presented with filled circles and triangles. The LAEs of triangles have a broad-band {V 
or i') magnitude fainter than that of the 2a level, and show the lower limits of the narrow- 
band excess color. The red and blue open circles denote spectroscopically-identified objects 
in the redshift range of LAEs and interlopers, respectively. We define the redshift ranges of 
LAEs as z = 3.09 - 3.18, 3.63 - 3.75, and 5.60 - 5.78 in the diagrams of NB503, NB570, 
and A^i?816, respectively (see Figure E] for these redshift ranges). The red and black squares 
mark AGNs from our multi-wavelength data (Table H]) with and without a spectroscopic 
redshift, respectively. For a display purpose, we place symbols at the narrow- v.s. broad- 
band color of 2.8 for objects with a color redder than 2.8. The green lines indicate 3a errors 
of the colors of narrow v.s. broad bands for a source with a color of V — NB503 = —0.10, 
V — NB 570 = 0.04, and i' — NB816 = 0.16, which correspond to the average color of all 
objects. The solid and dashed lines represent the 2a limit of a broad band, and the color cut 
of narrow-band excess. Note that the narrow-band magnitudes, i.e. A^-B503, NB570, and 
NB816, are total magnitudes, while the colors of ^ - NB503, V - NB570, and i' - NB816 
are defined with a 2"-diameter aperture. 
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Fig. 8. — Two color diagrams of continuum and narrow-band excess. Left panel: The 
diagram of B — V v.s. V — NB503 for our z = 3.1 LAEs. The black dots present colors of 
all the detected objects. The photometrically selected LAEs are plotted with filled circles 
and triangles. The LAEs of triangles have an off-band {V) magnitude fainter than 2a 
level, and show the lower limits of — A^i?503. The red and blue circles/triangles indicate 
spectroscopically identified objects in and out of the LAE redshift range, respectively. We 
take the LAE redshift range same as that of Figure [71 The red squares mark spectroscopically 
identified AGNs found in our multi- wavelength data (Table H]). For a display purpose, we 
place symbols at V — NB503 = 3.3 for the objects with a color of V — NB503 > 3.3. The 
objects with no continuum (V"-band) detection below a 2a level are plotted at a color of 
B — V = 0.6. The error bars at the upper-right corner indicate la errors for a typical LAE 
whose brightness is the median magnitudes of our photometrically selected LAEs. The solid 
line shows the color criteria of eq. ([T]) for sources with a detection in the broad band (V"). The 
colors are defined with a 2"-diameter aperture. The curves present tracks of model galaxies 
at different redshifts. Pink lines indica te model-LAE SEDs whic h are composite spectra of 
a 0.03 Gyr single burst model galaxy iBruzual fc Gharlot I l2003l and a Lyman a emission 
(i?PVi.cst(Lyct)=22 A); from the left to right, three different ampli tudes of IGM absorption 
are applied: 0.5tcs, Teg, and l.Srcff, where Teg is the iMadaul 's (Il995l ) original median opacity. 
The narrow-band excess in each of the peaks in the red lines indicates the Lym an a emission 
of LAEs at 2; = 3.1. Green lines show 6 templates of nearby starburst galaxies (IKinney et al. 



19961 ) up to z = 2, which are 6 classes of starburst galaxies with different dust extinction 
{E{B — V) = 0.0 — 0.7). The narrow-band excess peaks in the green line correspond to the 
emission lines of [Olll](z = 0.004), H/?(z = 0.03), or [Oil] (2: = 0.3). Black lines show colors 
of typical elliptical, spiral, and irregular galaxies (IColeman. Wu. &: Weedmanlll980l ) which 
are redshift ed fr o m z = to 2; = 2. Yellow star marks show 175 Galactic stars given by 
Gunn fc Strykerl (119831 ). Middle panel: Same as the left panel but for z = 3.7 LAEs. Right 
panel: Same as the left panel but for z = 5.7 LAEs. We do not show the error bar of a. R — i' 
color at the upper-right corner, because a typical LAE has no R-hand ffux above a la noise 
level. 
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Fig. 9. — Surface densities of objects detected in the NB503 data. The lower and upper 
points indicate surface densities of our z = 3.1 LAEs and all the objects detected in the 
narrow band, respectively. Black circles, hexagons, triangles, squares, and pentagons plot 
the surface densities of a ~ 0.2 deg field, SXDS-C, SXDS-N, SXDS-S, SXDS-E, and SXDS-W, 
respectively. We distinguish between the raw and completeness-corrected surface densities 
with the open and filled symbols, respectively. The red filled circles represent the surface 
density averaged over our ~ 1 deg^ survey field. The errors are given by Poisson statistics for 
black symbols, while the errors of red symbols are made by the geometric mean of Poisson 
errors and field-to-field variation of eq. (jlj). For the presentation purpose, we slightly shift 
all the black points along the abscissa. The exact magnitude of a black point is the same as a 
magnitude of neighboring red point. The vertical axis on the right side indicates the numbers 
of objects, i.e. iV/(0.5mag), identified in our ~ 1 deg^ survey area. The sohd, dotted, and 
dashed lines present the surface densities of LAEs reproduced by our simulations for the 
cases of a = —1.5, —1.0, and —2.0, respectively (see Section \5A\ for more details). 
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Fig. 10.— 



Same as Figure [HI but for NB570 data. 
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Fig. 11. — Same as Figure [9l but for A^i?816 data. The blue and cyan squares are the 
surface densities of z = 5. 7 LAE s det ected with the s ame NB816 filter but in different sky 
areas by lShimasaku et al.l (120061 ) and IHu et al.l (120041 1 . respectively. 



- 60 - 



d 




d 
<] 



0.0 




0.5 - 



1 .5 2.0 2.5 3.0 3.5 
V-NB503 



1.0 2.0 3.0 4.0 5.0 
V-NB570 



1.0 2.0 3.0 4.0 
i'-NB81G 



Fig. 12.— Color histograms of V" - A^5503 (left), V - A^S570 (middle) and i' - NB816 
(right) for the photometric sample of 2; = 3.1, 3.7, and 5.7 LAEs, respectively. The black- 
solid histogram presents the average of all the LAEs in the ~ 1 deg^ field. The shaded 
region shows the variations of histograms of five ~ 0.2 deg fields, SXDS-C, SXDS-N, SXDS- 
S, SXDS-E, and SXDS-W. The red histogram is the same as the black-solid histogram, 
but with a completeness correction. The error bars of the red histogram are given by the 
geometric mean of Poisson errors and field-to-field variation of eq. (jlj). The right-hand 
vertical axis on each panel indicates the numbers of LAEs for the ~ 1 deg^ histogram. We 
plot color distribution reproduced by our simulations of fixed a = —1.5 with dotted, dashed, 
solid, and dot-dashed lines, whose ctew values are presented in the legends. 
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Fig. 13. — Spectra of LAEs with AGN activities. We show four AGN-LAEs that have 
spectroscopic data (see Table S]). Dotted hnes with the legend indicate wavelengths of a 
typical emission from AGN. "©" marks the wavelengths where significant residuals of sky 
subtraction remain due to strong OH sky lines on a faint continuum. 
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Fig. 14. — The spectral energy distribution (SED) of LAEs with (an) X-ray, infrared, 
and/or radio counterpart (s). The legend shows object names and the corresponding symbols. 
Filled and open symbols represent an object with and without a spectroscopic identification, 
respectively . Solid and dotted c urves are th e temp la tes of radio-qui e t and radio-loud AGNs 
taken from IeIvIs et al.l fll994h : iTelfer et all Jioi); iRichards etaP toO^ (see text). The 
thick and thin curves indicate the template SEDs normalized with the detection limits of 
Lja and X-ray, respectively. The horizontal solid bars present our detection limits in X-ray, 
infrared, and radio bands. 
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Fig. 15. — Comparison of Lya fluxes measured with our images, /Lya(phot), and spectra, 
/Lya(spec), for our LAEs at 2; = 3.1 (left), 3.7 (middle), and 5.7 (right). The circles and 
diamonds indicate LAEs with FOCAS and VIMOS spectra, respectively. Since our FOCAS 
observations used a relatively narrow slit, we apply slit-loss corrections for the Lya fluxes of 
FOCAS. The error bars of circles and diamonds show the ranges of a 95% confidence level. 
The solid line at each panel presents the equality of /Lya(phot) and /Lya (spec). We plot 
dotted lines for the differences from the equality by a factor of 2. 
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Fig. 16. — Luminosity functions (LFs) of LAEs at z = 3.1. The filled circles are LPs 
obtained by the classical method. The error bars include the field-to-field variation (see eq. 
H]). The open circles, hexagons, triangles, squares, and pentagons represent the LFs of five 
~ 0.2 deg fields, SXDS-C, SXDS-N, SXDS-S, SXDS-E, and SXDS-W, respectively, which 
are also derived by the classical method. In order to avoid the overlaps of points, we slightly 
shift all the open symbols along the abscissa. The exact luminosity of an open symbol is the 
same as a magnitude of a neighboring filled circle. The solid line is the best-fit Schechter 
function with a = — 1.5 estimated by our simulations, but not the Schechter function fitted 
to the LFs of the classical method. 
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Fig. 17. — Same as Figure [T6l but for z = 3.7. 
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Fig. 18. — Same as Figure [T6l but for z = 5.7. 
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Fig. 19. — The best measurements of Lya luminosity functions (LFs) at 2; ~ 3 (top), 4 
(middle), and 6 (bottom), together with the previous measurements. The cyan, blue, and 
red filled circles are the LFs calculated with the classical method at z = 3.1, 3.7, and 5.7, 
respectively, that are labeled as O08 in the legend. The corresponding best-fit LFs from the 
simulations are plotted with colored solid lines. These filled circles and solid lines are the 
same as those presented in Figures [iGlfTSl Note that the z = 3.1 LFs (cyan hues) are overplot- 
ted on the middle and bottom panels for references. The other symbols and dotted-lines are 
the previous measurements and b est-fit Schechter functions w hose references are indicate d 
egend on each p a nel: G 07- 



in the 

KOO- 



Kudritzki et al. 



Muravama et al 



Gronwall et all (120071) t;05— Ivan Breukelen et all (120051 ) 



(l2000f ) 
J2007h . 506- 



C98- 



Tapken et all (120061). A 06— lAiiki et all (120061). HOA 



LFs of 



Hu 



et al 



Aiiki et al 



Hu et al 



Cowie fc 



iul JlQQsh . £'0 5— IPawsonl (2m\. M07- 



Shimasaku et 




WOQ- 



Westra et all (120061 ). T06— 
Hu et all (2004), if04r— recalculated 



(120041 ). A04— lAiiki et all (I2004J). M 4— iMalhotra fc R,hoadsl (120041 ). A03 



(120031) . and i?01— iRhoads fc Malhotral (120011 ). We recalculate the LFs of 



(I2OO4I ) data, and obtain the po i nts of HOAr whose values are close to those of 



the similar re-estimation by iTapken et al.l (120061 ). The arrows in the top and middle panes 
represent the luminosity range where LFs are dominated by LAEs with AGN activities (see 
Section H]6]) . For the reader's eye guide, we plot ticks of SFR obtained from eq. ([9]) on the 
upper abscissa axis. 
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Fig. 20. — Error ellipses of our Schechter parameters, L^ya and 0*. The left panel shows 
the ellipses for our observational data, while the right panel presents those for all LAEs with 
a positive emission {EW > 0) estimated by our simulations. Blue, green, and red contours 
represent z = 3.1, 3.7, and 5.7 LAE-LFs with the fixed slope of a = —1.5. Thick and thin 
lines indicate 1 and 2a confidence levels, respectively. 
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Fig. 21. — Inferred intrinsic Lya luminosity functions (LFs), together with the best-estimate 
apparent Lja LFs. Dashed- and sohd-hne curves are the intrinsic and apparent Lja LFs, 
respectively. Cyan, blue, and red curves denote for LFs at 2; = 3.1 (top panel), 3.7 (middle 
panel), and 5.7 (bottom panel), respectively. Cyan hues are repeatedly plotted at each panel 
for comparison. Light-gray and gray shades around cyan lines indicate 1 a errors of intrinsic 
and apparent Lya LFs aX z = 3.1. Areas of vertical and hatched lines are the same, but for 
z = 3.7 (blue) and 5.7 (red) in the middle and bottom panels, respectively. 
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Fig. 22. — The UV luminosity functions (LFs) of LAEs at z = 3.1 — 5.7 The cyan, blue, and 
red circles indicate the best-estimates of LFs for our z = 3.1 (top panel), 3.7 (middle panel), 
and 5.7 (bottom panel) LAEs. The solid lines represent the best-fit Schechter functions 
for the best estimates. The best-fit Schechter function oi z = 3.1 LAEs is presented on 
each panel for reference. The triangles are the lower limits of these LFs at each redshift 
that are derived with the aperture photometry of broad-band images (see text). For circles 
and triangles, we show the reliable and less-reliable measurements with filled and open 
symbols, which are obtained with UV magnitudes at > 5cr and 2 — 5cr levels, respectively. 
We only use these > 5a data (filled circles) for our Schechter- function fit (see text). The 
fill ed squares and st a rs on the b ottom panel are th e UV LFs of LAEs at z = 5.7 obtained 
by IShimasaku et al.l (120061 ) and IHu fc Cowid (120061 ) . On each panel, we also plot the UV 
LFs of dropout galaxies f or con iparison. The dash ed and dot-dashed lines are the z = 3 
LFs of Steidel et al. Jl999h and (jPaltani et al.ll2006f) . respectively. The dotted line indicates 



Steidel et all (Il999l ). but 0* of the LF is muhiplied by 1/10. Note that 



the z = 3 LF of 
the LF of z = 4 dropout galaxies are a. l most same as that o f z 



3 (ISteidel et all Il999 



Ouchi et al.l l2004al : iBeckwith et al.l l2006l : lYoshida et al.ll2006l ). On the bottorn pane l, we 
show th e UV LFs of dropout galaxies at ^ ~ 6 with crosses (IBouwens et al.l 120061 ) and 
asterisk (IShimasaku et al.ll2006l ). The gray region on the bottom panel indic ates the range 



of th e best-fit Schechter funct i ons for i-dropouts ob t ained by various stud ies (IBunker et al. 



2004J : lYan fc WindhorstI l2004l : iMalhotra et al.l l2005l : iBouwens et al.ll2006f ). which show the 
uncertainties of measurements for z = 6 dropout LF. For the reader's eye guide, we plot 
ticks of SFR obtained from eq. (|TOl) on the upper abscissa axis. 
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Fig. 23. — Histogram of the rest-frame Lya equivalent widths (EW) for our LAEs at z = 3.1 
(top), 3.7 (middle), and 5.7 (bottom). The black histograms and curves represent the best- 
estimate EWo and the EWo probability distribution estimated from errors of measurements, 
respectively, for photometrically selected LAEs with logL(Lya;) > 42.6 and EWq^*' > 70 — 80. 
The gray histograms and curves are the same, but for all the photometrically selected LAEs. 
The ticks of vertical axes indicate a fraction (left-hand side) and the number (right-hand 
side) in a bin size of AEW = 50A. The right-hand ticks correspond to the number of all 
the photometrically selected LAEs. The red histograms present the EWq distribution of our 
spectroscopically identified LAEs. On the bottom panel, the magenta histogram and black 
dashed line show the best-estimate EWo and the EWq probability distrib ution for z = 5.7 



LAE sample obtained by DEIMOS spectroscopy in the Subaru Deep Field ( IShimasaku et al. 
2006[ l. 
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Fig. 24. — Same as Figure [23| but for the comparison of histograms for apparent ("No 
IGMcorr") and intrinsic ("IGMcorr") EWq from z = 3.1 to 5.7. At each panel, the gray 
region indicates the apparent EWo histogram with uncertainties whose area corresponds 
to the allowed parameters of the black curves and lines (+error bars) of Figure [231 The 
cyan, blue, and red areas are the same as the gray regions, but for intrinsic EWq corrected 
for the IGM absorption, assuming the average Lya opacity (see text). For comparison, we 
repeatedly plot the 2; ~ 4 gray and blue areas with the black and blue meshes, respectively, 
on the top and bottom panels. 
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Fig. 25. — Rest-frame equivalent width {EWq) as a function of UV magnitude for 2; ~ 3 
(left), 4 (middle), and 6 (right) objects. Filled and open circles plot our spectroscopically 
identified LAEs with a UV magnitude brighter and fainter than our 3cr limits, respectively. 
The filled circles marked with an open square indicate LAEs with A GN activities. On th e 
left panel, we show EWq values of z ~ 3 dropo ut galaxies ob t ained by lShapley et al.l (120031 ). 
which are the average of the four subsam ples oflSteidel et al.l (120031) . On the right panel, the 

The stars and diamonds 



"(I2OO6I). 



filled squares represent z = 5.7 LAEs of iShinaasaku et al. 

present EWn of z ~ 6 dropout galaxies from IStanway et al.l (j2007l ) and the co mpilation of 



Ando et al.l (120061 ). The crosses plot EWq of z ~ 5 dropout galaxies presented in lAndo et al. 



( 120061 ). The dotted lines are loci of the constant Lya luminosity for 10^^, 5 x 10^^, 2 x lO''^, 
10^^, 5 X 10"^^, and 1 x 10*^^ erg s~^ from top left to the bottom. The shaded regions present the 
magni tude range with EW-large emitters (Muv ^ —21.5), which is claimed by lAndo et al. 
(l2006h . 
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Fig. 26. — UV-continuum color as a function of UV magnitude at 2 ~ 3 (left panel) and 
2; ~ 4 (right panel). The filled squares indicate a. i' — z' color of our LAEs at z = 3.1 and 
3.7 from our spectroscopic samples. The filled circles are the same, but for a. R — i' color. 
For confused LAEs whose photometry appears to be contaminated by neighboring objects 
on a broad-band image, we use plus signs and crosses instead of filled squares and circles, 
respectively. We mark LAEs with AGN activities wit h a large open squar e. The diamonds 



plot i' — z' colors of spectroscopically identified LBGs JVoshida et aPbooeh. The star marks 



on the left panel represent the average colors of LBGs obtained by lShapley et al.l (120031 ). We 
estimate i' — z' colors of Shapley et al.'s data with E(B — and the eq. [121 The vertical 
axis on the right-hand side shows the UV-slope index, /3, that corr esponds to i' — z' colors 
at 2; ~ 4. Note that this /? i s based on the origi nal definition (See ICalzettill200ll ). which is 
different from defined by lOuchi et al.l (l2004al ) . The difference of /3 values between R — i' 
and i' — z' colors is only < 7% in the color range of —0.3 — 0.3. The E[B — F) values 
corresponding to the UV colors are indicated with dotted lines. The gray curves represent 
the size of ±la error for a continuum-fiat (/j, = const) object. 
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Fig. 27. — Composite spectra of our LAEs at {z) = 3.13 (left), 3.68 (middle), and 5.68 
(right). The dotted lines indicate the wavelengths of interstellar absorption from star forma- 
tion, or high ionized emission lines from AGN activities. The plots of spectrum arc magnified 
in the inset boxes for wavelength ranges of Lja (left) and HcII (right). For the (z) = 5.68 
LAE, we only show the inset box of Lya, since our spectrum does not cover the wavelength 
of Hell emission for z — 5.7 objects. 



